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DEFINITIONS 
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decisions  attecting  major  programs,  (b)  address  Issues  ot  signiticant  concern  to  the 
Executive  Branch,  the  Congress  and/or  the  public,  or  (e)  address  issues  that  have 
significant  economic  implications.  IDA  Reports  are  reviewed  by  outside  panels  ot  experts 
to  ensure  their  high  quality  and  relevance  to  the  problems  studied,  and  they  are  released 
by  the  President  of  IDA. 

Group  Reports 

Group  Reports  record  the  findings  and  results  of  IDA  established  working  groups  and 
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that  they  meet  the  high  standards  expected  ot  retereed  papers  in  professional  journals  or 
formal  Agency  reports. 
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substantive  work  done  in  quick  reaction  studies,  (b)  to  record  the  proceedings  of 
conferences  and  meetings,  (c)  to  make  available  preliminary  and  tentative  results  of 
analyses,  (d)  to  record  data  developed  In  the  course  of  an  Investigation,  or  (e)  to  forward 
Information  that  is  essentially  unenalyzed  and  unevaluated.  The  review  of  IDA  Documents 
is  suited  to  their  content  and  intended  use. 
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PREFACE 


This  paper  was  prepared  under  the  Institute  for  Defense  Analyses’  Central  Research 
Program.  It  is  a  derivation  and  expansion  of  some  attrition  equations  incorporated  into 
IDA's  NAVMOD  model  of  naval  combat  (Reference  [3]). 

1  he  authors  are  grateful  to  Dr.  Lowell  Bruce  Anderson,  Dr.  Ronald  Enlow,  Dr. 
Arthur  Fries,  and  Dr.  Robert  Gilmore  for  their  review  of  this  paper. 
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ABSTRACT 


This  paper  examines  a  class  of  combat  processes  where  one  side’s  resources  are 
evenly  distributed  along  a  barrier  and  the  other  side  attempts  to  penetrate  that  barrier. 
Starting  from  a  set  of  straightforward  assumptions,  the  paper  derives  exact  or  approximate 
expressions  for  the  expected  numbers  of  resources  destroyed,  the  expected  number  of 
shots  fired  at  a  target,  and  in  some  cases,  the  probability  distribution  of  the  number  of  shots 
fired  at  a  target.  The  derivations  combine  and  extend  several  concepts  used  in  previous 
attrition  models,  including  having  resources  situated  on  a  barrier,  letting  each  attacker 
engage  several  targets,  and  the  use  of  the  Poisson  approximation  to  the  binomial  when 
there  are  several  types  of  targets.  Several  different  firing  rules  (attack  protocols)  are 
considered.  Some  of  the  equations  developed  here  have  been  used  in  the  Institute  for 
Defense  Analyses'  NAVMOD  naval  model  to  model  attrition  when  a  task  force  or  group  of 
ships  crosses  a  hostile  submarine  barrier. 
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INTRODUCTION 


This  paper  examines  a  class  of  combat  processes  where  one  side's  resources  are 
evenly  distributed  along  a  barrier  and  the  other  side  attempts  to  penetrate  that  barrier. 
Starting  from  a  set  of  straightforward  assumptions,  the  paper  derives  exact  or  approximate 
expressions  for  the  expected  numbers  of  resources  destroyed,  the  expected  number  of 
shots  fired  at  a  target,  and,  in  some  cases,  the  probability  distribution  of  the  number  of 
shots  fired  at  a  target.  The  derivations  combine  and  extend  several  concepts  used  in 
previous  attrition  models  [5,  6,  8,  9,  10].  Some  of  these  concepts  are  1)  having  resources 
situated  on  a  barrier,  2)  letting  each  attacker  engage  several  targets,  and  3)  the  use  of  the 
Poisson  approximation  to  the  binomial  when  there  are  several  types  of  targets.  Several 
different  attack  protocols  (firing  rules)  are  considered. 

Some  of  the  equations  developed  here  have  been  used  in  the  NAVMOD  naval 
model  [3]  to  model  attrition  when  a  task  force  or  group  of  ships  crosses  a  hostile  submarine 
barrier.  However,  the  barrier  could  be  composed  of  any  type  of  resource  (e.g.,  aircraft) 
spaced  evenly  along  a  barrier  line.  Accordingly,  in  the  rest  of  the  paper  we  shall  speak  of 
"barrier  elements"  and  "penetrators"  with  the  understanding  that  these  could  possibly  be 
ships,  submarines,  aircraft  or  any  appropriate  types  of  resources.  In  the  processes 
examined  here,  there  can  be  several  types  of  penetrators  but  only  one  type  of  barrier 
element.  (Strictly  speaking,  there  can  be  several  types  of  harrier  elements  present  as  long 
as  the  effectiveness  parameters  considered  are  the  same  for  all  the  types.) 

Chapter  I  considers  the  cases  where  the  barrier  elements  are  the  shooters  and  the 
penetrators  are  the  targets.  In  Chapter  II,  the  penetrators  are  the  shooters  and  the  barrier 
elements,  the  targets.  The  issue  of  combining  two  one-sided  attrition  processes  into  a 
process  where  attrition  to  both  sides  is  computed  is  not  addressed  in  this  paper.  Reference 
[1]  gives  one  method  for  doing  this;  Reference  [7]  gives  some  justifications  for  this 
method.  This  paper  is  structured  so  that  in  Chapter  I,  Section  G,  which  contains  results, 
can  be  read  immediately  after  Section  A,  the  statement  of  the  problem,  without  loss  of 
continuity.  Similarly,  in  Chapter  II,  Section  E  can  be  read  after  Section  A  without  loss  of 


1 


continuity.  Appendices  A  and  B  provide  proof  details  of  the  results  presented  in  Chapters  I 
and  II,  respectively. 

The  following  notation  is  used  throughout  the  paper.  If  A  and  B  are  events,  then 
P(A)  denotes  the  probability  of  the  event  A  and  AB  denotes  the  intersection  of  the  events  A 
and  B.  A  random  variable  is  denoted  by  a  symbol  with  a  tilde  over  it,  a  realized  value  of 
such  a  random  variable  by  the  corresponding  symbol  without  the  tilde.  Thus  the  random 
variable  h  might  assume  a  value  h;  its  expectation  is  E(h)  and  for  a  given  value  h,  one  can 
speak  of  the  probability  P(h=h).  The  symbol  "i"  should  be  read  as  a  lower  case  "ell."  For 
any  nonnegative  real  x,  LxJ  denotes  the  integer  part  of  x,  i.e.,  the  greatest  integer  less  than 
or  equal  to  x,  and  <x>  denotes  the  fractional  part  of  x,  i.e.,  <x>  =  x  -  LxJ. 
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I.  PENETRATORS  ARE  TARGETS, 
BARRIER  ELEMENTS  ARE  ATTACKERS 


A.  THE  PROBLEM 

Suppose  that  there  are  n  types  of  penetrators  (n  could  equal  1)  with  Rj  penetrators 
of  type  j,  for  j=l,...,n.  Let 

R  =  £  Rj 
j=i  J 


denote  the  total  number  of  penetrators.  There  is  one  type  of  barrier  element;  there  are  B 
barrier  elements  present.  We  consider  five  attrition  processes,  indexed  by  attack  protocol 
(a)  through  (e)  in  Assumption  5)  below,  that  proceed  according  to  the  following 
assumptions. 

1 )  Barrier  elements  are  spaced  evenly  along  a  barrier  of  length  L,  i.e., 
the  line  of  length  L  is  divided  into  B  equal  segments  and  a  barrier 
element  is  positioned  at  the  center  of  each  segment. 

2)  At  a  given  time,  each  penetrator  picks  a  crossing  point  x  along  the 
barrier  line,  according  to  a  uniform  distribution,  and  all  penetrators 
attempt  to  cross  the  barrier  simultaneously. 

3)  The  crossing  points  of  different  penetrators  are  mutually  independent 
random  variables. 

4)  If  the  crossing  point  of  a  given  penetrator  of  type  j  lies  within  distance 
wj  /2  of  a  given  barrier  element,  that  barrier  element  will  detect  the 
penetrator  with  probability  dj.  Otherwise,  that  barrier  element  will 
not  detect  the  penetrator.  Each  barrier  element  detects  the  penetrators 
that  come  within  wj  /2  of  it  independently  of  one  another. 

5)  Of  the  penetrators  it  has  detected,  each  barrier  element  chooses  certain 
ones  to  attack  with  certain  numbers  of  shots  (salvos),  according  to  one 
of  the  attack  protocols  described  below.  In  protocols  (b),  (c),  and  (d), 
the  parameter  H  is  that  number  of  salvos  available  per  barrier  element. 

(All  barrier  elements  use  the  same  attack  protocol.) 

(a)  A  barrier  element  fires  exactly  one  salvo  at  each  penetrator  it 
detects. 
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(b)  Of  the  penetrators  it  has  detected,  each  barrier  element  chooses  one 
according  to  a  uniform  distribution  and  fires  one  salvo  at  it.  The 
barrier  element  performs  this  process  exactly  H  times.  The  target 
choices  for  successive  firings  are  independent. 

(c)  If  a  barrier  element  detects  m  targets  (penetrators),  it  fires  [.H/mJ 
salvos  at  each  target,  then  chooses  according  to  a  uniform  distribution 

H-m[H/mJ  of  the  m  targets  and  fires  one  additional  salvo  at  each. 

(d)  If  the  number  of  targets,  m,  that  a  barrier  element  detects  does  not 
exceed  H,  the  shooter  fires  exactly  one  salvo  at  each  target.  If  m  > 

H,  the  shooter  chooses  according  to  a  uniform  distribution  H  of  the 
m  targets  and  fires  one  salvo  at  each. 

(e)  Of  the  targets  a  barrier  element  has  detected,  it  chooses  one  target 
according  to  a  uniform  distribution  and  fires  exactly  one  salvo  at  it. 
(This  is  the  special  case  of  protocols  (b),  (c),  and  (d)  when  H  =  1. 

We  consider  it  as  a  separate  protocol  because  the  formulas  for  the 
desired  quantities  are  considerably  simpler.) 

(6)  Conditional  on  detection  of  that  penetrator,  the  decisions  of  different  barrier 
elements  to  attack  a  given  penetrator  are  independent  events. 

(7)  A  shot  (salvo)  fired  at  a  penetrator  of  type  j  kills  it  with  probability  kj. 

The  effects  of  different  shots  (fired  at  the  same  penetrator)  are  independent. 


We  wish  to  compute  the  following  quantities: 

-  The  expected  number  of  penetrators  of  type  j  killed. 

(j  =  1,  •  •  .,n) 

-  The  expected  number  of  shots  (salvos)  fired  at  a  given  penetrator. 

-  The  probability  distribution  of  the  number  of  shots  fired  at  a  given 
penetrator. 

The  probabilistic  derivations  of  the  equations  assume  that  integer  numbers  of 
combatants  are  present.  Strictly  speaking,  therefore,  B  and  all  the  Rj  should  be  integers. 
In  the  formulas  which  have  a  closed  form,  however,  it  probably  does  little  harm  to  use 
nonintegral  values  for  B  and  the  Rj5  if  the  values  are  at  least  1.  (In  a  determinisitc  combat 

model,  the  iterative  use  of  attrition  equations  can  yield  noninteger  numbers  of  combatants, 
so  attrition  formulas  that  can  be  evaluated  for  noninteger  numbers  of  resources  are 
desirable.)  Where  a  formula  is  expressed  as  an  indicated  sum  with  a  limit  of  R  or  R-l, 
interpolation  between  the  results  with  LrJ  and  LrJ  +  1  is  a  reasonable  procedure.  At  least 
one  of  the  Rj  should  be  nonzero,  and  B  should  be  nonzero. 
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In  some  formulas,  the  number  of  possible  shots  per  shooter  H,  which  is  used  in 
attack  protocols  (b),  (c),  and  (d),  is  the  limit  of  an  indicated  sum,  and  thus  should  be  a 
positive  integer— but  reasonable  values  for  H  might  frequently  be  such.  (Throughout  this 
paper  the  terms  engagement,  shot,  and  salvo  will  be  used  synonymously.)  Reliability  of 
shots  can  be  considered  in  determining  the  values  of  the  kill  probability  kj.  The  barrier 
length  L  should  be  strictly  positive,  and  the  "detectability  width"  wj  should  be  less  than  or 
equal  to  L,  for  all  j.  (In  the  "homogeneous"  case,  i.e.,  n=l,  w,  d,  and  k  are  used  instead 
of  wj,  dj,  and  kj,  respectively.) 

Exact  or  approximate  formulas  for  the  three  main  quantities  of  interest  have  been 
developed  for  each  attack  protocol,  although  in  some  cases,  it  has  been  necessary  to  specify 
further  restrictions  on  the  parameter  values,  as  discussed  in  Section  D.  There  are  two 
reasons  why  some  of  the  formulas  are  only  approximations.  The  first  arises  from  the  use 
of  the  Poisson  approximation  to  the  binomial.  To  develop  algebraically  tractable  formulas, 
we  have  assumed  in  some  cases  that  certain  binomial  random  variables  are  approximately 
Poisson  distributed.  Section  E  explains  where  in  the  derivations  these  assumptions  are 
made. 

The  second  source  of  approximation  is  the  "ignoring  of  edge  effects."  This  paper 
derives  formulas  for  the  desired  quantities  in  terms  of  the  probabilities  of  certain  events  that 
involve  a  single  penetrator  and/or  a  single  barrier  element.  The  paper  posits  certain 
expressions  for  these  probabilities  that  are  reasonable,  and  are  exact  in  certain  cases,  but 
might  not  be  exact  if  a  barrier  element  is  located  near  the  edge  of  the  barrier,  or  if  the 
crossing  point  of  a  penetrator  is  near  the  edge  of  the  barrier  (where  "near”  depends  on  the 
values  of  certain  parameters).  These  expressions,  and  the  formulas  developed  from  them, 
thus  ignore  edge  effects.1  Section  F  discusses  some  of  these  expressions  and  the  reasons 
why  they  are  not  exact.  Other  such  expressions  are  explained  in  relevant  parts  of  the  paper. 

Without  loss  of  continuity,  the  reader  can  skip  to  Section  G,  which  presents  the 
formulas  that  have  been  developed  for  the  three  quantities  of  interest.  Each  theorem  in 
Section  G  states  if  edge  effects  have  been  ignored  in  its  derivation,  and  states  restrictions 
on  the  parameters  as  necessary.  The  formulas  are  stated  as  approximate  equalities  if  the 
Poisson  approximation  to  the  binomial  has  been  used  in  the  derivation. 


1  Exact  expressions  for  the  relevant  probabilities  would  have  to  take  the  edge  effects  into  account,  and 
might  be  more  complicated  than  the  expressions  used  here. 
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B.  DISCUSSION  OF  ASSUMPTIONS 


The  first  assumption,  the  geometric  placement  of  the  barrier  elements,  is  the  feature 
that  distinguishes  this  paper  from  certain  other  work  on  attrition  equations.  (Strictly 
speaking,  for  the  derivations  in  this  chapter  to  hold,  it  is  not  necessary  to  assume  that  a 
barrier  element  is  always  positioned  at  the  center  of  its  segment.  It  is  merely  necessary  that 
the  barrier  element  can  potentially  detect  and  attack  all  type-j  penetrators--and  only  those 
type-j  penetrators— that  transit  within  distance  wj/2  of  the  center  point.)  Another  barrier 
penetration  model  [2]  assumes  somewhat  more  independence  between  the  barrier  elements. 

The  second  and  third  assumptions  are  a  simplification  of  reality:  a  group  of 
penetrators  might  cross  a  barrier  in  some  coordinated  fashion,  or  sequentially.  In  the 
NAVMOD  model,  some  of  the  formulas  developed  here  are  used  to  model  barrier 
crossings,  but  provisions  are  made  for  the  penetrating  ships  to  cross  along  a  "corridor 
width"  that  may  be  less  than  the  full  physical  barrier  length  L.  The  details  appear  in  [3]. 
Note  that  the  uniform  distribution  of  the  penetrators’  crossing  points  implies  that  a 
penetrator  cannot  take  advantage  of  knowledge  of  the  barrier  elements'  positions.  In 
particular,  penetrators  cannot  discern  or  exploit  gaps  in  the  barrier  coverage. 

Throughout  this  chapter  we  will  use  the  following  terminology.  If  the  crossing 
point  x  of  a  particular  penetrator  of  type  j,  penetrator  p,  is  within  distance  wj/2  of  the 
location  of  barrier  element  P  (which  can  be  thought  of  as  the  midpoint  of  the  Pth  interval  of 
length  L/B  from  one  end  of  the  barrier)  we  will  use  the  phraseology  "p  is  detectable  by  P," 
or  ”p  is  vulnerable  to  P”,  or  "P  can  potentially  detect  p."  If  P  can  potentially  detect  p,  then 
P  detects  p  ("P  actually  detects  p,"  "p  is  detected  by  P")  with  probability  dj.  The 
"detectability  width"  wj  could  be  thought  of  as  a  sweepwidth  against  type-j  penetrators 
(especially  if  dj  =1).  It  is  possible,  if  the  wj  are  large,  that  the  detectability  areas  of 
different  barrier  elements  overlap,  i.e.,  that  a  penetrator  is  vulnerable  to  more  than  one 
barrier  element.  Frequently,  however,  it  is  reasonable  to  set  the  parameter  wj  to  L/B  for  all 
j.  The  interpretation  of  this  is  that  e.  ch  barrier  element  is  assigned  its  own  "barrier  cell," 
and  can  attack  only  those  penetrators  that  come  within  its  cell.  We  have  tacitly  assumed 
that  a  barrier  element  is  capable  of  attacking  any  penetrator  it  detects.  Therefore,  wj  should 
not  be  set  so  large  that  a  barrier  element  detects  penetrators  that  are  too  far  away  to  attack. 

The  second  part  of  assumption  4)  can  be  stated  more  precisely  as  follows.  Given 
any  specific  barrier  element  p  and  penetrators  (of  whatever  types)  pj,  p2,....,pjc  that  are 

detectable  by  p,  let  Dj  denote  the  event  that  p  actually  detects  penetrator  pj.  The  second 
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part  of  assumption  4)  states  that  the  Dj  form  a  set  of  mutually  independent  events.  This 
assumption  is  of  key  importance  in  the  derivations  of  the  formulas. 

We  note  that  in  assumption  5),  the  choices  of  target(s)  to  be  attacked  do  not  depend 
on  the  types  of  penetrators  detected  by  a  barrier  element,  but  only  on  the  total  number  of 
penetrators  detected.  In  general,  deriving  attrition  formulas  where  the  target  choice  is 
dependent  on  target  type  is  a  difficult  probabilistic  problem.  For  a  solution  of  one  special 
case  of  this,  see  Reference  [11]. 

If  the  detectability  widths  for  different  barrier  elements  (the  portion  of  the  barrier 
with  wj/2  on  either  side  of  a  barrier  element)  overlap  (this  could  happen  if  wj  or  B  is  large), 
it  is  possible  that  two  or  more  barrier  elements  can  detect  the  same  penetrator.  Assumption 
6)  merely  states  that  those  barrier  elements  do  not  coordinate  their  attacks.  This  is  not 
unreasonable  if  one  assumes  there  is  no  communication  between  barrier  elements  (e.g., 
barrier  submarines).  Assumption  7)  has  been  used  in  most  of  the  IDA  work  on  attrition 
equations  [5,9];  it  implies  that  if  a  penetrator  of  type  j  receives  exactly  h  shots,  it  will 
survive  with  probability  (l-kj)h.  This  is  a  realistic  assumption  in  the  case  where  a  fired 
shot  is  lethal  if  it  hits  the  target,  but  has  some  probability  (due,  say  to  unreliability)  that  it 
does  not  hit,  and  the  reliabilities  of  different  shots  are  independent.  On  the  other  hand,  it 
may  be  that  multiple  shots  hitting  a  target  have  a  cumulative  effect.  For  example,  there  is 
the  situation  where  one  hit  does  no  damage  at  all  to  a  target,  but  two  hits  destroy  it.  To 
model  a  situation  with  this  property,  it  is  often  desirable  to  know  the  probability 
distribution  of  the  number  of  shots  fired  at  a  target.  In  this  paper,  we  derive  this 
probability  distribution  for  some  cases.  Thus,  in  computing  attrition,  assumption  7)  could 
be  relaxed  and  a  function  giving  the  probability  a  target  is  destroyed  if  h  shots  are  fired  at  it 
could  be  used  instead.  NAVMOD  [3]  has  provisions  for  modeling  situations  where  hits 
have  a  cumulative  effect. 

We  note  finally  that  the  case  of  multiple  penetrator  types— n  >  1—  is  essentially 
intractable  if  the  exact  formulas  are  desired.  In  this  paper  we  derive  the  exact  formulas  for 
the  case  where  n  =  1  and  approximations  for  the  multiple  penetrator  type  case.  This  will  be 
discussed  more  in  Section  E  below. 

C.  RELATION  TO  PREVIOUS  WORK 

The  attrition  equations  developed  in  this  paper  are  pan  of  a  series  of  "binomial" 
attrition  models.  Such  models  can  be  considered  as  discrete  time  analogs  of  Lanchester 
attrition  equations.  Initial  numbers  of  searchers  and  targets  and  certain  effectiveness 
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parameters  are  input.  Attrition  occurs  over  an  interval  of  time,  according  to  specified 
assumptions  concerning  detection,  attack,  and  kill.  Formulas  for  the  expected  number  of 
targets  killed  (and,  in  some  cases,  other  quantities  of  interest)  are  probabilistically  derived 
from  the  assumptions.  The  initial  work  at  IDA  on  binomial  attrition  processes  appears  in 
Reference  [5].  Reference  [10]  derives  the  equations  for  many  such  processes;  Reference 
[9]  presents  results  for  most  of  the  important  processes.  Derivations  of  equations  for 
certain  specific  binomial  attrition  processes  appear  in  [6],  [11],  [12],  and  [13]. 

The  main  difference  between  this  paper  and  the  above  references  is  the  explicit 
positioning  of  the  resources  along  a  barrier  line.  The  previously  developed  binomial 
attrition  models  considered  uncoordinated  forces  of  attackers:  different  attackers  were 
mutually  independent.  Here,  because  of  the  barrier  spacing,  different  attackers  (barrier 
elements)  are  not  independent.  (Different  penetrators  are  independent.)  To  derive  expected 
kills  here,  we  condition  on  the  crossing  point  of  a  penetrator  and  then  still  must  ensure  that 
certain  random  variables  are  conditionally  independent  given  this  crossing  point.  This  line 
of  argument  was  unnecessary  in  previous  work.  We  point  out  that  attack  protocol  (a)  (each 
attacker  shoots  once  at  each  target  it  detects)  was  trivial  to  treat  when  attackers  were 
uncoordinated.  Introducing  the  barrier  structure  makes  protocol  (a)  more  interesting  to 
examine. 

This  paper  also  treats  in  an  integrated  fashion  the  issues  of  multiple  shots  per 
shooter  and  multiple  target  types.  Of  the  references  discussed  here,  only  Reference  [6] 
allowed  attackers  to  have  several  shots:  attack  protocol  (d)  was  used  and  only  one  type  of 
target  (and  one  type  of  shooter)  was  considered.  The  attrition  equation  developed  in  [6] 
(which  is  also  reported  in  [9])  does  not  have  a  closed  form:  it  is  expressed  in  terms  of 
cumulative  binomial  probabilities.  The  same  is  true  of  most  of  the  equations  developed  in 
the  current  paper  for  attack  protocols  (b),  (c),  and  (d),  where  an  attacker  has  more  than  one 
shot  available. 

Binomial  attrition  processes  involving  targets  (and  shooters)  of  more  than  one  type 
have  been  explored  much  more  fully.  If  detection  parameters  are  a  function  of  target  type, 
the  expected  number  of  targets  of  type  j  killed  can  be  rigorously  derived,  but  the  formula  is 
essentially  intractable  algebraically  [5].  To  circumvent  this  difficulty,  the  Poisson 
approximation  to  the  binomial  has  been  used  to  develop  approximate  attrition  equations  that 
have  simpler  forms  [8].  In  some  cases,  error  bounds  on  the  approximation  accuracy  of  the 
final  attrition  equation  have  been  developed  [8,  10].  (Kill  parameters  can,  without 
difficulty,  be  functions  of  both  shooter  type  and  target  type.  Reference  [12]  examines  the 
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process  where  detection  parameters  are  functions  of  shooter  type  but  not  target  type.  The 
resulting  attrition  equation,  first  derived  in  [5],  has  a  relatively  simple  closed  form.) 

Two  other  models  that  consider  combat  processes  similar  to  the  ones  treated  in  this 
paper  are  "An  Attrition  Model  for  Penetration  Processes”  [2],  and  the  attrition  routines  of 
the  "Barrier  Air  Defense  Model"  [4,  especially  pp.  E-43ff).  Reference  [2]  examines  a 
combat  process  in  which  penetrators  attempt  to  get  through  a  set  of  defenders,  that  can  be 
considered  as  forming  a  barrier.  Probability  distributions — not  just  expected  values — of  a 
wide  variety  of  quantities  are  computed.  The  assumptions  of  [2]  are  somewhat  different 
from  those  of  this  paper:  penetrators  attack  one  by  one,  in  succession;  the  attack  can  result 
in  attrition  to  both  penetrators  and  defenders,  and  each  defender  detects  a  penetrator 
independently  of  other  defenders.  Reference  [2]  admits  that  there  are  situations  where  this 
latter  assumption  is  not  realistic,  particularly  if  the  defenders  are  evenly  spaced  barrier 
elements. 

The  Barrier  Air  Defense  Model  considers  a  group  of  evenly  spaced  penetrating 
aircraft  opposed  by  defending  aircraft.  Several  combat  processes  are  considered;  in  one, 
the  defending  aircraft  are  also  evenly  spaced,  in  others,  the  defending  aircraft  are  positioned 
to  maximize  expected  attrition  to  the  penetrators.1  A  defender  can  detect  only  those 
penetrators  that  come  within  the  defender's  "launch  opportunity  zone";  detected  penetrators 
can  be  attacked  according  to  a  number  of  protocols  similar  to  those  in  Assumption  5  (of  this 
paper).  The  assumption  of  evenly  spaced  (rather  than  independent  random)  penetrators 
distinguishes  Reference  [4]  from  the  work  reported  here.  For  cases  where  penetrators  tend 
to  cluster  when  crossing  a  hostile  barrier,  the  assumption  of  independent  random 
penetrators  provides  a  tractable  way  of  examining  the  effects  of  such  clustering.  A 
clustered  crossing  pattern  is  not  the  same  as  independent  random  penetrators,  but  the  latter 
case  is  tractable  mathematically  and  is  in  some  sense  conservative:  no  coordination  or 
communication  between  different  penetrators  is  assumed.  The  results  can  thus  provide  a 
benchmark  for  evaluating  the  effectiveness  for  the  penetrating  side  of  any  particular 
crossing  pattern. 

Reference  [14]  examines  a  penetration  process  that  involves  a  single  transiting  ship 
crossing  a  minefield.  Each  mine  is  placed  along  a  barrier  line  according  to  a  uniform 
distribution;  the  locations  of  the  mines  are  independent  of  one  another.  The  crossing  point 


1  In  cases  where  the  barrier  must  be  formed  before  the  positions  of  the  penetrators  are  known,  even  spacing 
of  barrier  elements  along  the  barrier  line  would  seem  to  be  the  more  appropriate  assumption. 
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of  the  penetrator  is  uniformly  distributed  along  the  barrier  line  and  is  independent  of  the 
positions  of  the  mines.  The  probability  that  the  transitor  is  killed  is  derived,  and  the 
influence  of  edge  effects  is  thoroughly  examined  (cf.  Section  F  of  this  chapter). 

D.  AN  INDEPENDENCE  CONDITION 

We  have  been  able  to  derive  formulas  for  certain  of  the  desired  quantities  only  if  it 
can  be  guaranteed  that  certain  random  variables  are  independent  or  conditionally 
independent.  In  addition  to  the  basic  assumptions  1)  through  7),  certain  constraints  on  the 
parameters  are  sometimes  necessary.  To  explain  the  reasons  for  this,  let  us  examine  how 
one  might  go  about  deriving  the  formula  for  the  expected  number  of  penetrators  of  type  j 
killed. 

This  quantity  is  the  number  of  penetrators  of  type  j  multiplied  by  the  probability  that 
a  specific  penetrator  of  type  j  is  killed  by  the  barrier  elements  [12].  In  symbols, 

E(ARj)  =  Rj(l-P(S)) , 

where  the  random  variable  ARj  is  the  number  of  penetrators  of  type  j  killed  and  S  is  the 
event  that  a  given  penetrator  (call  it  penetrator  p)  of  type  j  survives  the  barrier.  We  thus 

wish  to  compute  P(S).  Let  the  specific  barrier  elements  be  indexed  by  b,  ranging  from  1 
through  B.  For  each  b,  let  Sp  be  the  event  that  penetrator  p  survives  barrier  element  b. 
The  intersection  of  the  events  S^,  is  S,  but  the  individual  S5  are  clearly  not  independent:  if 
penetrator  p  survives  barrier  element  (3,  the  chance  is  increased  that  p's  crossing  point  was 
not  within  the  detectability  area  of  [3,  thus  the  crossing  point  is  more  likely  to  be  in  the 
detectability  area  of  some  other  barrier  element  (call  it  y),  and  thus  the  probability  that  p 
has  survived  y  decreases.  The  interesting  question  is  whether  the  events  Sb  are 
conditionally  independent  given  the  crossing  point  x  of  penetrator  p. 

Consider  the  case  where  the  detectability  widths  of  the  barrier  elements  overlap 
(i.e.,  wj  >L/B),  and  where  the  attack  protocol  is  (b),  (c),  (d),  or  (e),  so  that  the  number  of 

shots  available  per  barrier  element  is  limited.  Suppose  that  penetrator  p's  given  crossing 
point  x  is  such  that  p  is  vulnerable  to  several  barrier  elements  (call  them  (3 1 ,  P2»~->Pk)- 
The  probability  that  barrier  element  Pi  fires  at  p  depends,  in  part,  on  the  number  of  other 
penetrators  barrier  element  Pi  detects.  In  any  realization  of  the  process,  the  crossing  points 
of  the  other  penetrators  are  fixed.  If  the  penetrator  survives  barrier  element  Pi,  this  might 
indicate  that  fewer  sh  ms  than  average  were  fired  at  p  by  Pi— which  might  have  happened 
because  Pi  detected  many  penetrators  and  had  to  spread  its  fixed  number  of  shots  thinly. 
This  would  imply  that  many  penetrators  crossed  within  the  detectability  area  of  Pi,  which 
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would  indicate  that  the  number  of  penetrators  crossing  within  the  detectability  area  of  some 
other  barrier  element,  say  $2>  was  less  than  average.  Since  p  is  vulnerable  to  p2> 
chance  of  p  receiving  more  shots  from  p2*  hence  being  killed,  would  be  increased.1 

The  above  argument  merely  indicates  why  it  is  not  obvious  that  the  events  Sb  are 
conditionally  independent  given  the  crossing  point  x  of  penetrator  p;  it  does  not  imply  that 
the  Sb  are  not  conditionally  independent.  Unfortunately,  numerical  examples  have  been 
developed  that  show  that  the  events  Sb  are  indeed  not  necessarily  conditionally  independent 
given  x.  In  this  case,  P(S)  must  be  found  by  integrating  over  the  appropriate  joint 
probabilities  of  the  specific  crossing  points  of  all  penetrators,  an  involved  calculation  that 
we  have  not  attempted  in  this  paper. 

Under  certain  conditions,  however,  the  Sb  do  form  a  set  of  conditionally 
independent  events.  Specifically,  with  p,  x,  and  Sb  defined  as  above,  the  following  two 
results  hold.  (Proofs  are  in  Appendix  A.) 

Lemma  1.1 :  If  attack  protocol  (a)  is  used  (each  barrier  element  fires  exactly  one  shot  at 

each  penetrator  it  detects)  then  the  Sb  are  conditionally  independent  given  x. 

Lemma  1.2:  If  penetrator  p  is  of  type  j  and  wj  <  L/B,  then  the  Sb  are  conditionally 

independent  given  x.  This  is  true  whether  or  not  wi  <  L/B  for  other  target 
types  i. 

One  reasonable  scenario  where  the  premises  of  Lemma  1.2  are  met  occurs  when  the 
barrier  line  (of  length  L)  is  partitioned  into  B  equal  nonoverlapping  "cells"  and  a  barrier 
element  is  allowed  to  attack  only  those  penetrators  it  detects  whose  crossing  points  lie 
within  its  cell.  Then,  for  all  j,  the  parameter  wj  is  the  minimum  of  the  "true  effective  sweep 
width”  and  L/B. 

If  protocol  (a)  is  used  or  if  wj  <  L/B,  we  will  use  the  terminology  "the 
independence  condition  holds."  The  independence  condition  is  used  in  the  derivations  of 
the  expected  number  of  penetrators  killed  and  the  probability  distribution  of  the  number  of 
shots  fired  at  a  penetrator.  The  formulas  for  the  expected  number  of  shots  fired  at  a 
penetrator,  however,  are  valid  whether  or  not  the  independence  condition  holds. 


1  This  argument  is  similiar  in  spirit,  though  not  in  detail,  to  the  explanation  of  why  targets  do  not  die 
independently  of  one  another,  given  in  Reference  [5],  p.  5. 
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E.  DERIVATION  OF  THE  DESIRED  QUANTITIES 


This  section  explains  the  main  arguments  used  in  the  derivation  of  the  three  desired 
quantities,  presenting  the  formulas  when  attack  protocol  (b)  is  used.  Proof  details  appear 
in  Appendix  A.  The  homogeneous  (n  =  1)  and  heterogeneous  (n  >  1)  cases  are  considered 
separately,  where  appropriate.  In  the  course  of  the  derivation,  we  point  out  certain  issues 
of  approximation,  which  are  discussed  more  fully  in  Section  F. 

We  start  the  derivation  by  considering  one  specific  penetrator,  penetrator  p,  which 
is  assumed  to  be  of  type  j,  and  one  specific  barrier  element  (3  that  can  potentially  detect  p. 
We  then  use  the  results  of  this  situation  to  derive  the  overall  quantities  of  interest. 

1.  One  Penetrator,  One  Barrier  Element 

We  assume  that  the  crossing  point  x  of  penetrator  p  is  such  that  barrier  element  P 
can  potentially  detect  p-i.e.,  x  is  within  wj/2  of  the  location  of  P  (which  is  at  the  midpoint 
of  the  Pth  interval  of  length  L/B  along  the  barrier  line).  We  define  the  following  events  and 
random  variables: 

D-p  actually  detects  p 

Av— The  number  of  penetrators  other  than  penetrator  p  (regardless  of  type),  that  p 
detects,  is  exactly  i.  (Defined  for  i=0,l,...,R-l.) 

DAx— The  intersection  of  events  D  and  Ax,  i.e.,  P  detects  p  and  exactly  i  other 
penetrators. 

K— random  variable  indicating  the  number  of  shots  P  fires  at  p. 

F^-The  event  that  K=h,  i.e.,  P  fires  exactly  h  shots  at  p.  (Defined  for 
h=0, 1  ,...,H.) 

K-p  kills  p. 

We  compute  P(K),  E(h),  and  the  probabilities  P(Fh).  These  correspond  to  our 
three  desired  quantities,  but  in  the  case  where  only  one  barrier  element  and  one  penetrator 
are  present.  Later,  we  integrate  these  quantities  into  overall  results.  Before  we  derive 
P(K),  E(h),  and  the  P(Fh)  we  derive  their  "component  parts"  P(Fh  I  DAt),  P(K  |  Fh),  and 
P(DAV). 

Because  the  choice  of  targets  p  attacks  depends  only  on  the  number,  not  the  types, 
of  penetrators  detected,  the  conditional  probabilities  P(FhlDAj.)  are  functions  of  h  and  t 
only.  If  attack  protocol  (b)  is  used,  from  the  statement  of  Assumption  5b)  it  is  clear  that,  if 
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P  detects  p  and  i  other  targets,  i.e.,  i+l  targets  in  total,  the  number  of  shots  target  p 
receives  is  binomially  distributed  with  parameters  H  and  l/(i+l).  In  symbols. 


P(FhIDA)  = 


i+l 


(  V** 


) 


i+l 
V  J 


This  is  defined  for  h=0,...,H  and  i=0,...,R-l,  and  in  fact  is  sensible  mathematically  for 
any  nonn^gative  i.  (Formulas  for  P(FhlDAL)  for  the  other  attack  protocols  appear  and  are 
derived  in  Appendix  A  [Lemma  1.5].) 


Let  a  bar  over  any  event  denote  its  complement,  so  K  is  the  event  that  P  does  not 
kill  p.  From  Assumption  7),  it  is  clear  that 

P(KlFh)  =  (l-k/, 
so 

P(KIFh)  =  1  -  (1— kj)h,  for  every  h. 

The  probabilities  P(DAi)  are  somewhat  more  interesting.  Since  penetrator  p  is 
assumed  to  be  detectable  by  p,  P(D)  is  simply  dj.  Consider  some  other  penetrator,  p', 
which  is  of  type  i.  In  order  to  be  detected  by  P,  penetrator  p'  must  first  be  detectable  by 
P— that  is,  its  crossing  point  must  be  within  w[/2  on  either  side  of  the  location  of  P- 
Ignoring  edge  effects,  the  probability  of  this  is,  by  Assumption  2),  2(wj/2)/L,  or  wi/L. 
Given  that  p'  is  detectable  by  P,  p'  will  actually  be  detected  by  p  with  probability  di.  Thus 

P(penetrator  p',  of  type  i,  is  detected  by  p)  =  diwi/L. 

The  potential  edge  effects  problems  will  be  discussed  in  Section  F,  below.  By  Assumption 
3),  crossing  points  of  different  penetrators  are  independent,  and  by  the  second  part  of 
Assumption  4),  penetrators  detectable  by  barrier  element  P  are  actually  detected  by  P 
independently  of  one  another.  Thus,  if  pi,  P2>— ,Pm  are  different  penetrators,  the  events 

(P  actually  detects  penetrator  pk) 

for  k=l,...,m  form  a  mutually  independent  set.  By  this  same  reasoning,  the  event  D  (P 
detects  p)  is  independent  of  the  event  P  detects  pk,  for  pk  *p-  There  are  R  -  1 
penetrators  other  than  penetrator  p.  Therefore,  if  there  is  only  one  type  of  penetrator,  the 
number  of  penetrators  other  than  penetrator  p  that  barrier  element  P  detects  is  binomially 
distributed  with  parameters  R-l  and  wd/L.  In  symbols. 
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f. 


and 


P(A )  = 

i 


P(DA )  =  dP(A )  =  d 


R-l 

i 


'wdV  ' 


L 

V  J 


wd 


'  R-ll 

wd^ 

l 

wd^ 

l  1  J 

l  LJ 

l  -TJ 

\  R-l-i 


\  R-l-i 


for  i=0,...,R-l. 


Now  suppose  that  n  >  1.  (The  following  argument  is  identical  to  that  given  in 
Reference  [9].)  Recall  that  the  specific  penetrator  p  under  consideration  is  of  type  j.  Not 
including  penetrator  p,  there  are  R j  -  1  penetrators  of  type  j  and  Ri  penetrators  of  type  i,  for 
all  i=l,...,n,  i*j.  Let  the  random  variable  I-  denote  the  number  of  type-i  penetrators,  other 
than  penetrator  p,  detected  by  p.  It  is  clear  from  the  preceding  arguments  that: 

-  For  i*j,  Tj  is  distributed  binomially  with  parameters  Rj  and  wj  di  /L. 

-  i.  is  distributed  binomially  with  parameters  Rj  -1  and  wj  dj  /L. 

-  ij ,  are  a  mutually  independent  set  of  random  variables. 

The  sum 


i=l 


is  the  total  number  of  penetrators  other  than  penetrator  p  that  P  detects.  The  distribution  of 
T  is  a  convolution  of  binomial  distributions  with  different  parameters,  and  the  resulting 
formula  is  algebraically  intractable  [5].  We  therefore  use  the  Poisson  approximation  to  the 
binomial  and  consider  each  T  as  being  approximately  Poisson  distributed  with  mean 
RiWjdj/L,  for  i  *  j.  Random  variable  Ij  is  approximately  Poisson  with  mean  (Rj-l)wjdj/L. 
Since  the  sum  of  independent  Poisson  random  variables  is  Poisson,  T  is  approximately 
distributed  Poisson  with  mean 

n 

=  (Rj-DWjdj/L  +  XRjW^/L. 

»*j 

Then  the  probability  P(A0  =  P(t=i),  i.e.,  the  probability  that  barrier  element  p  detects 
exactly  t  targets  other  than  penetrator  p  is  approximately 


M2 


and 


od  e 


-°5 


A! 


P(DA  ) »  d=  *hl  1 

1  J  i! 

We  now  return  to  the  derivation  of  P(K),  E(h) ,  and  P(Fh).  The  events  Fo 
through  Fh  form  a  partition,  thus 

H  H 

P(K)  =  £P(KlFh)P(^)  =£P(KlFh)P(Fh), 

h=0  h=l 

since  if  (3  fires  zero  shots  at  p,  K  will  not  occur.  Similarly 

H 

E(h)  =  hP(Fh) . 

h=l 

The  probability  that  penetrator  p  receives  exactly  h  shots  from  barrier  element  (5  depends  on 

the  total  number  of  targets  (3  has  detected.  The  events  D,  DAo,  DAi,...,DAr_i  form  a 
partition,  so 

R-l 

P(Fh)  =  P(FhlD)P(D)  +  £  P(FhlDA)P(DA  )  . 

1=0  11  t  i 

We  have  assumed  that  no  shot  fired  by  a  barrier  element  can  reach  a  penetrator  not  detected 
by  that  barrier  element.  Thus  P(FqID)  =  1  and  P(FhlD)  =  0forallh>l.  Therefore 


H  R-l 

P(K)  =  £  £  P(KIFh)  P(FhlDA)P(DA  ) 

h=i  i=0  11 

H  R-l 

E(K)=£  £  hP(Fh  IDAi)  P(DAV)  . 

h=l  i=0 


The  appropriate  formulas  derived  earlier  can  be  substituted  for  the  indicated 
probabilities,  and  the  algebra  performed.  The  following  points  should  be  noted.  First,  it  is 

generally  easier  to  compute  P(K)  and  E(h)  by  interchanging  the  order  of  summation  and 
summing  on  h  First.  The  formulas  for  P(Fh)  are,  in  general,  complicated  algebraically. 
Second,  many  of  the  formulas  do  not  reduce  to  a  closed  form  and  must  be  left  as  indicated 
summations.  Third,  in  the  heterogeneous  case  (n>l)  the  approximate  (Poisson)  formulas 
for  P(DAt)  will  be  used.  To  compensate  for  the  Poisson  approximation  to  the  binomial. 


1-13 


indicated  sums  on  i  will  be  taken  over  all  nonnegative  values  of  i,  not  just  up  to  R-l. 
(Note  that  the  formulas  for  both  P(FhlDAt)  and  the  approximation  to  P(DAi)  are  well 
defined  for  all  nonnegative  integer  t,  and  Fubini's  theorem  justifies  interchanging  the  order 
of  summation  on  h  and  t.) 

In  attack  protocol  (b),  most  of  the  formulas  will  not  reduce  to  closed  forms.  Let  us 
introduce  the  following  alternative  notation  for  P(DAt).  In  the  homogeneous  case  (one 
type  of  penetrator),  we  denote  P(DAL)  by  where 


If  n  >  1,  we  first  define 

a 


and  then  the  Poisson  approximation  to  P(DAi)  by  \j/jl,  where 

=  dje^aj/i!  i  =  0,l„ 


w  d 

=  (R,-l)— ~ 


y 

,ti  •  L 
i*j 


This  assumes  that  the  specific  penetrator  p  under  consideration  is  of  type  j.  (A  further 
approximation,  which  is  somewhat  easier  computationally,  is  to  replace  aj  by 

-  Y  D  d‘W' 

«  =  A  Ri-f— 

1=1  L 


in  each  Reference  [8]  does  this;  here  we  will  leave  the  formulas  in  terms  of  the  aj.) 


Using  this  notation,  after  algebraic  simplification,  the  formulas  for  the  quantities 


P(K),  E(h),  and  P(Fh),  where  attack  protocol  (b)  is  used,  are  as  follows.  In  the  homogen¬ 
eous  case 


P(K)  = 

E(h)  = 


d  - 

LH 

wR 


f  \ 

i+l-k 


H 


V 


i+l 


) 


wdN 

R~ 

1  - 

1  u 

v  ^ 

J’ 
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R-l 

P(F0)  =  1-d  +  X  “ 

i=0 


and 


P(Fh)  =  £  co 


i=0 


fH] 

M 

ii 

(  \ 

t 

Is 

i+i 

V  J 

t+i 

vl+1y 


H-h 


In  the  heterogeneous  (n  >  1)  case,  assuming  that  the  specific  penetrator  is  of  type  j,  the 
following  approximate  formulas  obtain: 


P(K)  =  dj  -  X  V 

i=()  Ji 


x+1 


dH  -a 

E(h)  =  —  (1-e  J)  , 

a 

J 


P(F())  -  1  -  d,  +  ^  V. 


(  \H 

i 


vl+1y 


and 


P(Fh)  -  £  Vji 

i=0 


^  H 

vhy 


vl+1y 


f  \H-h 
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V  J 


h=I . H. 


2.  Combined  Results  Considering  All  Barrier  Elements 

The  previous  subsection  considered  one  penetrator  p,  of  type  j,  and  one  barrier 
element  p  that  could  potentially  detect  p.  Using  the  quantities  P(K),  E(h),  and  the  P(Fh) 
derived  in  Subsection  E.l,  we  now  derive  the  overall  quantities  of  interest.  We  again 
consider  one  penetrator  p,  assumed  to  be  of  type  j.  The  individual  barrier  elements  are 
indexed  by  b,  b=l,...,B.  We  use  the  following  notation,  some  of  which  was  introduced  in 
Section  D. 

Sb--event  that  barrier  element  b  does  not  kill  penetrator  p  (i.e.,  p  survives  b). 

g 

S-o  S. ,  i.e.,  the  event  that  p  survives  all  the  barrier  elements. 
b=i  b 

y  -random  variable  indicating  the  number  of  shots  Fired  at  p  by  barrier  element  b. 
b 
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B 


y~- 


b=l 


i.e.,  the  total  number  of  shots  fired  at  p 


x-point  on  the  barrier  at  which  penetrator  p  crosses,  uniformly  distributed  on 

[0,L]. 

Cx  -event  that  x  =  x,  defined  for  x  €  [0,L]. 

Nj(x)-number  of  barrier  elements  that  can  potentially  detect  p,  if  p  is  of  type  j 
and  the  crossing  point  of  p  is  x.  If  n=l  we  will  use  the  notation  N(x). 

Again,  the  quantities  we  wish  to  derive  are: 

-  The  expected  number  of  penetrators  of  type  j  killed. 

-  The  expected  number  of  shots  fired  at  a  given  penetrator. 

-  The  probability  distribution  of  the  number  of  shots  fired  at  a  given  penetrator. 

In  terms  of  the  above  notation,  these  desired  quantities  are  given  by  Rj(l-P(S)), 
E(y  )  and  P(y=y)  (for  appropriate  y),  respectively.  We  derive  each  quantity  in  turn. 

To  find  P(S)  we  first  condition  on  the  crossing  point  x.  We  shall  speak  of 
probabilities  conditional  on  Cx,  even  though  P(CX)  for  zero  for  every  specific  x.  The 

crossing  point  x  is  distributed  uniformly  on  [0,L],  thus 

L 

P(S)  =  J  P(SICx)j-dx. 

0 

If  the  independence  condition  holds,  as  described  in  Section  D,  then  the  events  Sb  are 

conditionally  independent  given  x,  so  that 

B 

PCS  ICX)  =  II  P(SbICx),  x  e  [0,L], 
b=l 

If  x  is  not  within  the  detectability  area  of  barrier  element  b,  the  penetrator  will  certainly 
survive  barrier  element  b,  so  P(SblCx)  =  1;  if  x  is  within  the  detectability  area  of  b,  then 
by  the  results  of  the  previous  section,  P(SblCx)  =  1  -  P(K),  regardless  of  the  specific 
barrier  element  b  involved.  Thus 

B  N  (x) 

II  P(SJC  )  =  (l-P(K))  J 
b=l  b  x 
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The  geometric  spacing  of  the  barrier  elements  gives  rise  to  the  following  formula  for 
Nj(x). 

WjB 

Lemma  1.3:  Let  I.  denote  the  integer  part  and  f-  the  fractional  part  of  —  •  tin 
J  J  L 

w.B  w.B  w.B 

symbols  I.  =  — j- —  and  f.  =  - ~ —  .)  Then,  ignoring  edge  effects ,  Nj(x)  equals 

J  [  L  J  J  L  L  L  J  J 

Ij  on  a  set  of  measure  (l-fj)L,  and  equals  lj  +  1  on  a  set  of  measure  fjL.  Alternatively,  if 


the  random  variable  Nj(x)  is  the  number  of  barrier  elements  that  can  potentially  detect 
penetrator  p,  where  x  is  the  uniformly  distributed  crossing  point  of  penetrator  p,  then 


N  (x)  = 

J  1  !j+1 


Proof:  See  Appendix  A.  Section  F,  below,  discusses  the  effect  of  edge  effects  on 
Nj(x). 

Combining  the  above  results,  we  obtain 

I  F+l 

PS)  =  (l-fjMl-PCK))3  +f}(\-P(K))}  , 

which  simplifies  to 

I. 

P(S)  =  (1-P(K))J  (l-fjP(K)). 

We  note  that  P(K)  is  itself  a  function  of  the  type  j  of  penetrator  p  .  The  expected  number  of 
penetrators  of  type  j  killed  is 

E(A~Rj  )  =  Rj  (l-P(S)) 

=  Rj[l-(l-P(K))Ij(l-fjP(K))]. 

The  computation  of  E(y)  proceeds  in  a  similar  manner.  We  again  consider  a 


specific  penetrator  p,  of  type  j.  Conditioning  on  the  crossing  point  x  of  p,  we  obtain 

fL  1 

E(y)  -  J  E(y  I C)  y  dx- 
0  x  L 
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For  each  x. 


E(y!Cx)  =  X  E(yblcx). 

b=l 

If  x  is  such  that  barrier  element  b  cannot  potentially  detect  p,  then  b  with  certainty  fires  zero 
shots  at  p.  If  b  can  potentially  detect  p,  then  by  the  results  of  Subsection  1,  the  expected 

number  of  shots  b  fires  at  p  is  the  quantity  denoted  by  E(h) .  Therefore 

E(y!Cx)  =  Nj(x)E(K). 

This  is  true  whether  or  not  the  independence  condition  holds.  Integrating  over  x  and 
recalling  the  formula  for  Nj(x),  we  obtain 


Noting  that  I.  +  f.  1= 


E(y)  =  (1-fj)  IjE(h)  +  fj(Ij+l)  E(K). 
v/B 

— s—  t  we  can  simplify  this  to 


W:  B 

E(y)  =  ~~  E(h) 


HThis  formula  makes  sense  in  that  the  average  number  of  barrier  elements  that  can  potentially 
detect  p  is  the  average  value  of  Nj(x),  or 

w.B 

(1-f.)  I.  +  f.  (I.  +  1)  =  I.  +  f.  =  . 

J  J  J  J  J  J  L 

and  each  barrier  element  that  can  potentially  detect  p  fires  an  average  of  E(K)  shots  at  p. 

To  compute  the  probability  distribution  of  y,  we  must  assume  that  the  premises  of 


the  independence  condition  hold,  i.e.,  attack  protocol  (a)  is  used  or  wj  <  L/B.  If  protocol 
(a)  is  used,  each  barrier  element  fires  exactly  one  shot  at  each  penetrator  it  detects.  We 
have  made  the  tacit  assumption  that  detections  of  a  given  penetrator  by  different  barrier 
elements  are  conditionally  independent  given  the  crossing  point  x  of  that  penetrator.  (They 
are  clearly  not  independent  overall.)  Under  protocol  (a),  then,  if  penetrator  p  crosses  at  x, 
ihe  number  of  shots  fired  at  p  is  binomially  distributed  with  parameters  Nj(x)  and  dj.  That 


is 


P(y=yiCx) 


,  y  , 


d'.  (1-d.) 
j  j 


Nj(x)— y 
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for  y  =  0,l,...,Nj(x).  Integrating  over  x,  we  obtain,  after  extensive  algebra 


V 


y  lry 

dI 


l-fJ  + 


fjd-dj)  (Ij+1) 

(Ij  +  1  -  y) 


P(y=y}  =<  fdV1 


J  J 


0 


y=o . ij 

y=ij+i 


for  all  other  y. 


If  Wj  <  L/B,  then  by  Lemma  1.3,  Nj(x)  is  either  0  or  1,  and  penetrator  p  is  never 
vulnerable  to  more  than  one  barrier  element.  If  Nj(x)  =  0,  then  no  barrier  element  can 
potentially  detect  p,  and  y  is  zero  with  certainty.  If  Nj(x)  =  1 ,  then  one  barrier  element, 
call  it  (3,  can  potentially  detect  p,  and  y  has  the  same  distribution  as  h,  which  was  derived 
in  Section  E.l.  If  wj  <  L/B,  fj  =wj  B/L,  and  Nj(x)  is  zero  with  probability  1  -  fj.  Thus 


P(y=0)  =  (1-f.)  •  1  +  f  P(F0) 
w.B 

=  1  -  -~t~  (1-P(F0) ) 

and 

P(y=y)  =  (l-fj)  •  0  +  fj? (Fy) 

w.B 

=  -j-  P(Fy)  ,  y=i . H. 


It  can  easily  be  verified  that  these  formulas  are  also  valid  when  wj  =  L/B,  i.e.,  Ij  =  1  and 
fj  =0. 

F.  EDGE  EFFECTS 

Edge  effects  can  affect  the  accuracy  of  the  formulas  for  the  three  quantities  of 
interest  (the  expected  number  of  penetrators  killed,the  expected  number  of  shots  fired  at  a 
penetrator,  and  the  probability  distribution  of  shots  fired  at  a  penetrator)  if  some  barrier 
elements  are  located  less  than  wj/2  away  from  the  end  of  the  barrier,  for  some  penetrator 
types  j.  In  this  case,  the  "detectability  width"— the  interval  of  length  wj  centered  on  the 
location  of  the  barrier  element— overlaps  the  area  beyond  the  end  of  the  barrier,  where  no 


penetrators  cross.  In  terms  of  the  symbols  defined  in  the  preceding  section,  this  has  two 
major  consequences: 

1.  The  true  value  of  Nj(x)~the  number  of  barrier  elements  that  can  potentially 
detect  a  penetrator  of  type  j  crossing  at  point  x--is  lower  than  the  formula  given 
previously.  See  the  proof  of  Lemma  1.3  in  Appendix  A  for  more  discussion  of 

Nj(x). 

2.  For  barrier  elements  b  near  (i.e.,  less  than  wj/2  away  from)  the  edge  of  the 
barrier,  the  probability  that  a  given  penetrator  is  vulnerable  to  barrier  element  b 
is  not  wj/L  but  equals  1/L  times  (wj/2  plus  the  distance  of  b  from  the  barrier 

edge),  a  smaller  quantity.  This  changes  the  probabilities  P(AO  and  P(DAi) 

used  in  computing  P(K),  E(h),  and  P(h=h). 


It  is  not  immediately  obvious  how  these  two  effects  affect  the  accuracy  of  the 

formulas  for  the  overall  quantities  E(AR.),  E(y),  and  P(y=y).  However,  the  problem  of 

edge  effects  may  not  be  exceedingly  important.  If  barrier  elements  are  repositioned  so  that 
no  element's  detectability  area  extends  beyond  the  ends  of  the  barrier,  the  values  of  both 
Nj(x)  and  P(AV)  will  increase,  and  the  errors  in  these  two  quantities  might  compensate  for 
each  other.  And  if  the  independence  condition  holds,  edge  effects  are  often  automatically 
eliminated,  as  described  below. 


By  Assumption  1),  if  there  are  B  barrier  elements,  the  barrier  line  is  divided  into  B 
equal  segments  and  a  barrier  element  is  positioned  at  the  midpoint  of  each  segment.  For  a 
barrier  of  length  L,  if  one  endpoint  is  labeled  with  the  coordinate  value  0,  then  the  position 
of  barrier  element  b  is 

L  L 
xb  “  B^b  “  7 

for  b=l,  .  .  .,  B. 

In  particular,  barrier  elements  1  and  B  are  positioned  at  L/2B  and  L(1-1/2B), 
respectively.  Recall  from  Section  D  that  the  events  Sb  must  be  conditionally  independent 
given  the  crossing  point  x  of  a  penetrator  if  the  formulas  for  and  E(ARj)  and  P(y=y)  are  to 
be  valid.  (Here,  "valid"  means  that  a  formula  is  correct,  possibly  subject  to  Poisson 
approximation  or  edge  effects.)  This  independence  condition  will  hold  if  either 


(a)  Attack  protocol  (a)  is  used  (a  barrier  element  fires  exactly  one  shot  at 
each  penetrator  it  detects) 

(b)  wj<L/B  for  the  particular  penetrator  type  j  considered. 
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or  both.  As  stated  earlier,  frequently  wj  will  be  less  than  or  equal  to  L/B  for  all  j, 

specifically  in  the  case  where  each  barrier  element  is  not  allowed  to  attack  penetrators 
outside  the  barrier  element's  "cell"  of  length  L/b.  Therefore  the  formulas  for  E(A~Rj ),  E(y) 

and  P(y=y)  are  valid  for  all  j.  But  also  in  this  case,  no  barrier  element  is  less  than 
wj/2  away  from  the  barrier  edge,  for  any  j,  and  thus  edge  effects  are  not  a  problem  either. 

The  issue  of  what  happens  when  wj<L/B  for  some  but  not  all  j  is  more 
complicated.  Consider  penetrators  of  type  j  where  wj<L/B.  The  independence  condition 

does  hold,  and,  except  for  edge  effects,  the  formulas  for  E(AlL) ,  E(y)  and  P(y=y)  are 

valid.  We  conjecture,  however,  that  they  are  not  subject  to  edge  effects  if  and  only  if  attack 
protocol  (a)  is  used.  The  problem  is  that  (in  the  derivation  presented  in  this  paper)  for  other 
attack  protocols  we  must  condition  on  the  number  of  other  penetrators  a  barrier  element 
detects-and  the  appropriate  probabilities  of  detection  will  be  subject  to  edge  effects. 

Less  holds  true  for  penetrators  of  those  types  j  for  which  wj>L/B.  If  attack 
protocol  (a)  is  used  the  independence  condition  holds,  and  the  formulas  for  E(ARj),  E(y), 
and  P(y=y)  are  valid  but  subject  to  edge  effects.  If  an  attack  protocol  other  than  (a)  is  used, 
the  formulas  for  E(Altj)  and  P(y=y)  are  not  valid,  and  the  formula  for  E(y)  is  subject  to 
edge  effects. 


Edge  effects  will  become  more  and  more  of  a  problem  the  greater  the  ratios 
wj/(L/B)  are.  For  the  distance  of  barrier  element  b  from  the  "left-hand"  barrier  edge  (point 

0)  is 


and  this  will  be  less  than  wj/2  away  from  this  edge  if 


2b  <  1  + 


L/B 


Thus  the  greater  the  ratio  wj/(L/B),  the  more  barrier  elements  b  satisfy  this  inequality  and 
thus  have  their  detectability  widths  overlap  the  area  beyond  the  left-hand  barrier  edge. 
Similar  arguments  apply  to  the  right-hand  barrier  edge  (point  L). 
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G.  THE  RESULTS 


This  section  presents  the  formulas  for  the  three  main  quantities  of  interest— the 
expected  number  of  penetrators  killed,  the  expected  number  of  shots  fired  at  a  given 
penetrator,  and  the  probability  distribution  of  the  number  of  shots  fired  at  a  given 
penetrator.  Results  are  given  for  the  five  attack  protocols  in  both  homogeneous  (n=l)  and 
heterogeneous  (n>l)  cases.  Results  for  which  we  have  not  been  able  to  find  a  closed  form 
are  expressed  in  terms  of  cumulative  binomial  or  Poisson  probabilities  when  possible.  The 
methods  of  derivation  of  these  formulas  have  been  presented  in  Section  E,  above;  proof 
details  appear  in  Appendix  A. 

Note  that  if  H  is  set  equal  to  1  in  the  formulas  for  attack  protocols  (b),  (c),  or  (d) 
(and  the  appropriate  algebraic  simplifications  are  performed),  the  corresponding  formula 
for  attack  protocol  (e)  results.  Also,  if  H  is  extremely  large,  the  formulas  for  protocol  (d) 
reduce  to  the  corresponding  formulas  for  protocol  (a).  The  descriptions  of  the  attack 
protocols  imply  that  the  above  results  should  hold;  that  they  do  is  an  independent  check  on 
the  accuracy  of  the  formulas. 

1.  Notation 

This  section  briefly  reviews  some  of  the  notation  used  in  the  formulas.  The 
quantities  for  which  formulas  are  presented  are  as  follows: 

E(ARj  )-expected  number  of  penetrators  of  type  j  killed.  (In  the  homogeneous 
case,  simply  E(£R).) 

E(y)-expected  number  of  shots  fired  at  a  specific  penetrator.  (In  the  hetero¬ 
geneous  case,  the  penetrator  is  assumed  to  be  of  some  specific  type  j.) 

P(y=y)--the  probability  that  exactly  y  shots  are  fired  at  a  specific  penetrator 
(of  type  j). 

The  following  parameters  are  defined  in  the  assumptions  in  Section  A:  L,  wj,  dj, 
H,  kj.  In  the  homogeneous  case,  w,  d,  and  k  are  used  instead  of  wj,  dj,  and  kj, 
respectively.  There  are  n  types  of  penetrators,  Rj  penetrators  of  type  j,  R  penetrators  in 
total  (R  =  Z  Rj),  and  B  barrier  elements. 

The  probabilities  that  a  given  barrier  element  detects  a  given  detectable  penetrator 
and  exactly  t  other  penetrators  (the  P(DAt))  appear  in  some  of  the  expressions  that  do  not 
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have  closed  forms.  For  conciseness,  we  will  use  the  notation  0)^  and  \j/ji  for  them.  In  the 
homogeneous  case  we  set 

f  R-l) 

0)  =  d  , 

i  V  ) 

We  will  often  denote  the  quantity  (wd/L)  by  p.  In  the  heterogeneous  case  we  first  define, 
for  each  j=l,...n, 


wdY  f  wd^  1-1 

1=0 . R-'- 


d-w 

a)  (Rj~1)_irL 


i=l  L 

i*j 


and  then  set 


_CJ  1 

V.  =d.  i=0,  1,2,...  . 

J  i  J  t ! 


The  following  notation  will  be  used  for  binomial  and  Poisson  probabilities.  (In  this 
paragraph  only,  n  is  any  nonnegative  integer  and  p  is  any  real  number  between  zero  and 
one.)  We  define 


b(m;n,p) 


£(m;n,p)  = 


n  m/i  d1-"1  n 

m  P  0~P)  ,  m=0,...,n 

0  for  m<-l 

m 

Z  b(i;n,  p)  for  m=0,...,  n 

t=0 

1  for  m>n 


and 


P(m ;  (i)  = 


-a  m  ,  , 

e  p  /m! 


f  0  for  m<-l 

P(m  ;  p)  =  j  Zp(m;  p)  for  m=0,  1,  2 

l  t=0 
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Finally  we  will  use  the  quantities  Ij  and  fj  defined  (in  the  heterogeneous  case)  for 
j=l,...,n  by 


In  the  homogeneous  case  we  use  the  analogous  quantities  I  and  f. 


2.  The  Expected  Number  of  Penetrators  Killed 
a.  Homogeneous  Case 

Theorem  I.la.  Suppose  that  n=l  and  attack  protocol  (a)  is  used.  Then, 
ignoring  edge  effects. 


E(AR)  =  R 


l-(l-kd)  (1-fkd) 


Theorems  1.1  h.  c.  d.  and  e.  Suppose  that  n=l  and  w<L/B.  E(AR)  is 


given  by 


E(AR)  =  R. 


ji 


1-(1-P(K»  (l-fP(K»j , 

where  P(K)  is  defined  for  the  specific  attack  protocol  used,  as  indicated  below.  (Note  that 
1=0  or  1,  and  if  1=1,  f=0.) 


-  Attack  protocol  (b) 


R-l 


?(K)  =  d  -  y  co.  i 

t+1 


H 


1=0 


-  Attack  protocol  (c)  Define  U  =  min(H-l,  R-l),  and  let 


U 

A  -XtMl-k) 

i=0 


_H 

i+l 


o^k) . 
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Then  if  H>R,  P(K)  =  d-A,  and  if  H<R-1 


P(K)  =  ^ 


+  dB(H-l;R-l,  p)  -  -™£kB(H;R,  p), 

wR 


where  p= 


Let  p=wd/L.  Then  if  H>R,  P(K)  =  kd,  and  if  H<R-1, 


P(K)  =  kdB(H-l;R-l,  p)  +  ^  (1-B(H;R,  p  ) ). 


-  Attack  protocol  (e) 


P(K)  l  >  ]- 


b.  Heterogeneous  Case 


Suppose  that  n>l  and  that  attack  protocol  (a)  is  used.  Then, 


ignoring  edge  effects, 


E(AR.)  =  R.  [l-(l-k.d  )  J  (1-f.k.d.)  ] . 

J  J  J  J  J  J  J 

Theorems  I.2b.  c.  d.  and  e.  Suppose  that  n>l  and  wj<L/B  for  all 

i=l,...,n.  Then  E(ARj)  for  j=l,...,n,  is  approximately  given  by 

I. 

E(ATL)  =  [1— (1— P(K))J  (l-f.P(K))] 

where  P(K)  is  defined  for  the  specific  attack  protocol  used,  as  indicated  below.  (Note  that 


1=0  or  l,and  if  I.  =1,  f.=0.) 
J  J  J 


Attack  protocol  (b) 


OO  f  V 

P(K)  -  d  -  £  M 


Attack  protocol  (c)  Define 


H-l  JL 

a  =  i  \t/(i-k)Lt+iJ(i^JL>k). 

1=0  J L  J  l+l  J 
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Then  if  aj>0 


and  if  aj=0 


P(K)  - 


Hk.d. 
J  J 
a. 

J 


Hk.d. 

+  d.P(H-l;  a.)  -  -±*- 


P(K)»d.[l-(l-k.)H] 


P(H;  a.) 


(It  is  assumed  that  H  is  a  positive  integer.) 


-  Attack  protocol  (di 

k.d.H 

k.d.P(H-l;  a.)  +-J-L.  (1-P(H;  a.))  ifa.>0 
P(K)  « <  J  J  J  “j  J  J 

k.d.  if  a.=0 
.  J  J  J 


Attack  protocol  (el 


P(K)  H 


k.d.  -a. 

-^(1-e  J)  if  a  >0 
j  1 


k.d.  if  a.=0. 
L  jj  j 


3.  The  Expected  Number  of  Shots  Fired  at  a  Given  Penetrator 

In  this  subsection,  we  consider  a  specific  penetrator,  which  we  call  penetrator  p.  In 
the  heterogeneous  case,  j  denotes  p's  type.  The  random  variable  y  is  the  expected  number 
of  shots  fired  at  penetrator  p  (summed  over  all  barrier  elements).  Here  we  present 
formulas  for  E(y);  in  subsection  4,  we  present  formulas  for  P(y=y),  for  appropriate  y. 


a.  Homogeneous  Case 

Theorems  I.3a.  b.  c.  d.  and  e.  Suppose  that  n=l.  Then,  ignoring  edge 
effects,  if  attack  protocol  (a),  (b),  (c),  (d),  or  (e)  is  used,  as  indicated  below,E(y)  is  given 
by  the  following  formulas: 

-  Attack  protocol  (a) 
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-  Attack  protocol  (o') 


E(y)  =  M[,_(i_M)R| 

E(y)  = 


-  Attack  protocol  (d~) 


If  H>R,  E(y)  =  as  in  (a). 

_  wBrrlTO/T4_1. 


If  H<R-1,  E(y)  =  ^°[dB(H-l;R-l,p)  +  — (1-B(H;  R,  p))], 
where  p=wd/L 


E(y)  =  |(1-<1-^)R] 

Note  that  attack  protocols  (b)  and  (c)  have  the  same  formula.  (See  Appendix  A  for  more  on 
this.)  If  H=l,  it  can  be  verified  that  the  formula  in  (d)  (and  of  course  those  in  (b)  and  (c)) 
reduces  to  that  in  (e). 

b.  Heterogeneous  Case 


ignoring  edge  effects, 


i.  Suppose  that  n>l  and  attack  protocol  (a)  is  used.  Then, 


E(y)  = 


Theorems  I.4b.  c.  d.  and  e.  Suppose  that  n>l.  Then,  ignoring  edge 
effects,  if  attack  protocol  (b),  (c),  (d),  or  (e)  is  used,  as  indicated  below,  E(y)  is 
approximately  given  by  the  following  formulas. 


E(y)  = 


w-BHd:  -a: 

- -  (1-e  0  if  a:>0 

T  ~  J 


Vj  BHd/L 


if  aj  =0. 
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f  w-BHdj 

- (1-e  J)  ifa->0 

E(y)  J  L“i 


Wj  BHdj/L 


if  a}  =0. 


as  in  (b). 

-  Attack  protocol  (d) 


w.B  Hd. 

-4—  [d.P(H-l;  a.)  +  — ^{1-P(H;  a.))  ]  ifa>0 


L  1  J 


w.Bd./L 
J  J 


Attack  protocol 


E(y)  - 


WjBdj  -<*j 
J 

w.Bd/L 
J  J 


if  a.=0. 
J 


if  a>0 
J 

if  a -0 
J 


4.  The  Probability  Distribution  of  the  Number  of  Shots  Fired  at  a  Given 
Penetrator 


a.  Attack  Protocol  (a) 

Theorem  I.6a.  Suppose  that  n>l  and  attack  protocol  (a)  is  used.  Then, 
ignoring  edge  effects  and  assuming  that  the  given  penetrator  is  of  type  j, 

rrn  ij-y  I"  fjd-dj)  (ij+i)" 

y  d/d-d/  >-fi+ - ,  +,.y  .  y=0 . Ij 

p(y=y)  =  1  ^  J  L  J  J 


y=ij  +1- 


Theorem  1. 5a.  Suppose  that  n=l  and  attack  protocol  (a)  is  used.  Then, 
subject  to  edge  effects,  the  formulas  for  the  probabilities  P(y=y)  are  as  given  in  Theorem 
1. 6a,  but  substituting  d,  I,  and  f  for  dj,  Ij,  and  fj,  respectively. 
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b.  Attack  Protocols  (b),  (c),  (d),  (e) 

Before  presenting  the  results  for  P(y=y)  for  these  cases,  we  relate  the  P(y=y)  to  the 
probabilities  P(Fh),  where  Fh  is  the  event  that  a  specific  penetrator  (penetrator  p)  receives 
exactly  h  shots  from  a  specific  barrier  element  (3  that  can  potentially  detect  p.  The  random 
variables 

hb=number  of  shots  fired  at  penetrator  p  by  barrier  element  b 

for  b=l,...,  B  are  in  general  not  independent,  and  not  even  conditionally  independent  given 
the  crossing  point  x  of  p.  To  avoid  problems  with  this,  we  assume  that  wj<L/B  for  all 
penetrator  types  i=l,.„,  n,  so  that  penetrator  p  is  detectable  by  at  most  one  barrier  element. 
The  following  lemma,  proved  in  Appendix  A,  is  valid  in  both  the  homogeneous  and 
heterogeneous  cases;  in  the  former  case,  merely  suppress  the  subscripts  i  and  j.  As  before, 
j  denotes  the  type  of  the  specific  penetrator  p  under  consideration.  For  protocols  (a)  and 
(e),  consider  H  as  1. 

Lemma  1.4:  Suppose  wi<L/B  for  all  i.  Then,  for  all  attack  protocols, 

-If  w}<L/B 

P(y=0)  =  l-fj  +  fjP(Fo) 

P(y=y)  =  fj  P(Fy ).  y=i . h 

-  If  Wj  =  l/b 

p(7=y)  =  P(Fy )  y=0 . H. 

Recall  that  Ij  and  fj  are  defined  as  the  integer  and  fractional  parts,  respectively,  of 
wjB/L.  If  wj<L/B,  then  lj=0  and  fj=wjB/L.  If  wj=L/B,  then  formally  Ij=l  and  fj=0,  but 
letting  Ij  =0  and  fj=l  in  the  first  set  of  formulas  yields  the  second  set  of  formulas. 
Accordingly,  for  the  remainder  of  this  section  only,  we  redefine,  for  j=l,...,n, 

[  wj  B 

I  -j—  if  w.<IVB 

f,  -  L 

y  1  if  Wj  =L/B. 

In  the  homogeneous  case,  we  define  f  by  the  above  formulas,  suppressing  the  subscript  j. 
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Utilizing  derived  expressions  for  the  P(Fy),  we  can  state  the  following  theorems  for 
the  homogeneous  and  heterogeneous  cases,  respectively. 

Theorems  T.5h.  c.  d.  and  e.  Suppose  n=l  and  w<L/B.  Then  the  random 
variable  y  has  the  probability  distribution 

P(y=0)  =  1  -  f  +  fP(F0) 

P(y=y)  =  fP(Fy).  y=l,...,  H, 

where  the  P(Fy)  are  defined  separately  for  attack  protocols  (b),  (c),  (d),  and  (e),  as 


indicated  below: 


Attack  protocol  (b). 


P(Fo)  = 


R-t  fnl  f  i  Y  ( 
P«V>  -  2>.  y  TTT  T+T 


y=i,...,H. 


Let  [xl+  denote  the  positive  part  of  x,  i.e., 

+  /x  if  x>0 

W  =  \0  if  x<0. 


and  define  p=wd/L.  Then 


P(F0)  =  l-d+d(l-B(H-l;  R-l,  p)) - —  (1-B(H;  R,  p  ) ). 


(Note  that  if  H>R,  this  formula  becomes  1— d,  which  is  appropriate  considering  the 
meaning  of  H,  R,  and  Fq.) 


P(F,)  =  2d(B(H-l;  R-l,p)-B( 


H 

yJ-l;R-l, 


-^(B(H;  R,  p) -B(  lyl;  R,  p  ) ) 


+  _Sr(1"B(H;R’p)) 


and  if  H>2,  then  for  y=2,...,H 


H 

H 

P(Fy)  =  (y+l)dlB( 

.  y. 

-1;  R-L  p)  -  B( 

.y+1J 

HL 

wR 


; 

Hi 

H 

[B(j 

Ly. 

;  R,  P)  -  B( 

Ly+1J 

;  R,  P  )  ]+ 


+ 


HL 

wR 


H 

H 

[  B( 

y-!_ 

;  R,  p  )  -  B( 

.  y. 

H 

H 

(y-l)  d  IB( 

Ly— ij 

-1  ;R-1,p)  -  B( 

.y. 

-  Attack  protocol  (d).  Let  p=wd/L.  Then 


where 


P(F0)  =  1  -P(Fi) 


P(Fj)  =  d  B(H-1;  R-l,  p)  + 


HL 

wR 


(1-B(H;  R,  p)). 


P(Fy)  =  0,  y>2. 


Note  that  if  H>R,  P(Fi)  becomes  d,  which  yields  the  same  result  as  protocol  (a)  and  is 
consistent  with  the  definitions  of  H,  R,  and  Fi. 


-  Attack  protocol  (e).  Again,  let  p=wd/L.  Then 

P(F0)  =  l-P(Fi), 


where 

p(f,>  =  i-<i-<i-p)R>. 

P(Fy)  =  0,  y>2. 

It  can  be  verified  that  if  H=l,  the  formulas  for  P(Fo)  and  P(F  1 )  in  attack  protocols  (b),  (c), 
and  (d)  reduce  to  those  in  protocol  (e). 
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Theorems  I.6b.  c.  d.  and  e.  Suppose  that  n>l  and  wj<L/B  for  i=l . n. 

If  the  specific  penetrator  p  under  consideration  is  of  type  j,  then  the  random  variable  y  has 


the  distribution 


P(y=0)  =  1-fj+fj  P(Fq) 

P(y=y)  =  fjP(Fy) »  y=l,...,H, 

where  the  P(Fy)  are  approximately  given  by  the  expressions  below,  defined  separately  for 
attack  protocols  (b),  (c),  (d),  and  (e). 

-  Attack  protocol  <b) 


P(F0)  .  M,+  | 


Vi 


P<Fy)-I 

t=0 


i+l 


-)y( 


,»-y 

) 

i+l  h 


y-1 . H. 


-  Attack  protocol  (c).  Let  [x]+  denote  the  positive  part  of  x,  as  explained  in  Theorem  I.5c. 
Then  if  aj  >0 


p(F0) 

P(Fj) 


Hd. 

Jf 


-  1-dj-Kl:  [1— IP(H— 1  ;a.)  ]  -  — ^l-P(H;a.)  ] 

J  J  J  j  J 


H 


«  2d.[P(H-l;a.)  -  P(|y 

Hdj  r  J  H| 

-  -^[P(H;a  )  -  P(jy;a)] 

J  J  J 


-Ua.)  ] 


Hd. 

+  -oT-H-PCH;^)] 
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and  if  H>2,  then  for  y=2,...,H 


H 

H 

P(Fy)  ~  (y+l)d.[P( 

_y. 

— l;otj)  -  P( 

Ly+1J 

Hd;  I  _Hj  *  + 

-  [P(|  "tn  ;  a. )-  IP  ^5  a.)] 


a. 

j 

Hd. 


a. 

j 


[P( 


y 

H 


y-1 


-  (y — l)  dj  [P( 


-«;o  )  - 


If  aj=0,  then  P(Fo)=l-dj,  P(Fn)=dj,  and  P(Fy)=0  for  all  other  values  of  y. 
-  Attack  protocol  (d'l 

P(Fq)  =  l-PCFi), 

where 


P(Fi) 


Hd 

d-  P(H-1;  a.)  +  — ri(l— P(H;  a-) )  ifa,>0 

j  j  a.  J  j 

j 

dj 


P(Fy)  =  0,  y>2. 

-  Attack  protocol  (el 

P(F0>  -  \ 


-a 


P(Fi)  -  < 


if  a-  =0. 


1 — -(l-e  J) 
ai 

N 

if  ctj  >0 

if  aj=0 

dj  -«■ 

—  (l-e  *) 

<aj 

if  otj  >0 

A 

ifaj=0 

II.  BARRIER  ELEMENTS  ARE  TARGETS, 
PENETRATORS  ARE  ATTACKERS 


We  now  examine  the  attrition  to  barrier  elements  during  the  barrier  crossing 
process.  The  notation  is  similar  (and  sometimes  identical)  to  that  of  Chapter  I.  We  still 
speak  of  barrier  elements  and  penetrators.  In  this  chapter,  however,  the  terms  "barrier 
element"  and  "target"  will  be  used  synonymously,  and  the  terms  "penetrator",  "attacker", 
and  "shooter"  will  be  used  synonymously.  The  statement  of  the  problem  in  Section  A, 
below,  is  almost  identical  to  the  corresponding  statement  in  Chapter  I,  but  the  penetrators 
and  barrier  elements  have  reversed  shooter  and  target  roles.  Some  of  the  implications  of 
this  are  discussed  in  Section  B,  below. 

A.  THE  PROBLEM 

Suppose  that  there  are  n  types  of  penetrators  (n  could  equal  1)  with  Rj  penetrators 
of  type  j.  Let 

r  -  iRj 

j=i 


denote  the  total  number  of  penetrators.  There  is  one  type  of  barrier  element;  there  are  B 
barrier  elements  present.  We  consider  five  attrition  processes,  indexed  by  attack  protocols 
(a)  through  (e)  in  Assumption  5)  below,  that  proceed  according  to  the  following 
assumptions. 

1)  Barrier  elements  are  spaced  evenly  along  a  barrier  of  length  L,  i.e.,  the  line  of  length  L 
is  divided  into  B  equal  segments  and  a  barrier  element  is  positioned  at  the  center  of  each 
segment. 

2)  At  a  given  time,  each  penetrator  picks  a  crossing  point  x  somewhere  along  the  barrier 
line,  according  to  a  uniform  distribution,  and  all  penetrators  attempt  to  cross  the  barrier 
simultaneously. 

3)  The  crossing  points  of  different  penetrators  are  mutually  independent  random  variables. 

4)  If  the  crossing  point  of  a  given  penetrator  of  type  j  lies  within  distance  wj/2  of  a  given 
barrier  element,  that  penetrator  will  detect  the  barrier  element  with  probability  dj. 
Otherwise,  that  penetrator  will  not  detect  that  barrier  element.  Each  penetrator  detects 
the  barrier  elements  that  lie  within  wj/2  of  the  penetrator's  crossing  point  independently 
of  one  another. 


5)  Of  the  barrier  elements  it  has  detected,  each  penetrator  chooses  certain  ones  to  attack 
with  certain  numbers  of  shots  (salvos),  according  to  one  of  the  attack  protocols 
described  below.  In  protocols  (b),  (c),  and  (d),  the  parameter  Hj  is  the  number  of 
salvos  available  per  penetrator  of  type  j.  (All  penetrators  use  the  same  protocol.) 

(a)  A  penetrator  fires  exacdy  one  salvo  at  each  barrier  element  it  detects. 

(b)  Of  the  barrier  elements  it  has  detected,  each  penetrator  chooses  one  according  to  a 
uniform  distribution  and  fires  one  salvo  at  it.  The  penetrator  performs  this  process 
exacdy  Hj  times.  The  target  choices  for  successive  firings  are  independent. 


(c) 


If  a  penetrator  detects  m  targets  (barrier  elements)  it  fires 


salvos  at  each 


target,  then  chooses  according  to  a  uniform  distribution  H;  - 1  Hl  I  rn  of  the 

J  LmJ 

targets  and  fires  one  additional  salvo  at  each. 


m 


(d)  If  the  number  of  targets,  m,  that  a  penetrator  detects  does  not  exceed  Hj,  the 
shooter  fires  exactly  one  salvo  at  each  target.  If  m  >  Hj,  the  shooter  picks 
(randomly  and  uniformly)  Hj  of  the  m  targets  and  fires  one  salvo  at  each. 

(e)  Of  the  targets  it  has  detected,  a  penetrator  chooses  one  target  according  to  a  uniform 
distribution  and  fires  exactly  one  salvo  at  it.  (This  is  the  special  case  of  protocols 
(b),  (c),  and  (d)  when  Hj  =  1.) 

6)  Conditional  on  detection  of  that  barrier  element,  the  decisions  of  different  penetrators  to 
attack  a  given  barrier  element  are  independent. 


7)  A  shot  (salvo)  fired  by  a  penetrator  of  type  j  at  a  barrier  element  kills  the  barrier  element 
with  probability  kj.  The  effects  of  different  shots  (fired  at  the  same  barrier  element)  are 
independent. 


We  wish  to  compute  the  following  quantities: 


-  The  expected  number  of  barrier  elements  killed. 

-  The  expected  number  of  shots  (salvos)  fired  at  a  given  barrier  element. 

-  The  probability  distribution  of  the  number  of  shots  fired  at  a  given  barrier  element. 

Exact  or  approximate  formulas  for  the  first  two  of  these  quantities  have  been 
developed  for  each  attack  protocol.  General  formulas  for  the  probability  distribution  of  the 
number  of  shots  fired  at  a  given  barrier  element  have  been  derived  only  for  attack  protocols 
(a),  (d),  and  (e).  As  in  Chapter  I,  we  have  sometimes  ignored  edge  effects  when  deriving 
these  formulas.  Throughout  this  chapter  and  Appendix  B,  we  indicate  the  derivation  steps 
and  resultant  expressions  where  edge  effects  have  been  ignored.  Section  D  discusses  some 
of  the  consequences  of  ignoring  edge  effects. 
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Without  loss  of  continuity,  the  reader  can  skip  to  Section  E,  which  presents  the 
formulas  that  have  been  developed  for  the  three  quantities  of  interest. 

B.  COMPARISON  WITH  CHAPTER  I 

The  notation  for  the  problem  statement  has  been  chosen  to  reflect  the  comple¬ 
mentary  nature  of  the  combat  processes  examined  in  the  two  chapters.  The  parameters  n, 

Rj,  R,  B,  and  L  denote  exactly  the  same  quantities  in  Chapter  II  as  in  Chapter  I,  and  the 
same  restrictions  noted  in  Chapter  I,  Section  B  apply  to  them.  The  parameters  wj,  dj,  Hj, 
and  kj  denote  in  Chapter  II  quantities  similar  to  the  corresponding  ones  of  Chapter  I,  but 
they  apply  to  penetrators  rather  than  barrier  elements. 

In  this  paper,  penetrators  operate  (mutually)  independently  of  each  other,  but  barrier 
elements  are  positioned  in  a  fixed  pattern  (evenly  spaced  along  a  line).  There  can  be 
several  types  of  penetrators  but  just  one  type  of  barrier  element.  In  the  current  chapter,  the 
penetrators  are  the  shooters,  so  different  shooters  operate  independently.  This  makes  it 
easier  than  in  Chapter  I  to  derive  the  expected  number  of  targets  killed  and  the  expected 
number  of  shots  fired  at  a  given  target.  No  special  independence  condition  is  necessary. 
Since  there  is  only  one  target  type  there  is  no  need  to  use  the  Poisson  approximation  to  the 
binomial  to  obtain  tractable  formulas  for  the  number  of  targets  a  shooter  detects.  (This 
approximation  is  used  in  computing  the  probability  distribution  of  the  number  of  shots 
fired  at  a  target  in  the  case  of  heterogenous  shooters.)  It  is  relatively  easy  to  deal  with 
multiple  shooter  types.  Even  so,  the  similarity  of  many  of  the  final  formulas  to  the 
corresponding  ones  in  Chapter  I  is  striking. 

We  use  the  words  "detectability,"  "vulnerability,”  et  al.,  in  a  manner  similar  to 
Chapter  I.  A  penetrator's  "detectability  width"  or  "detectability  area"  is  the  line  of  length  wj 
(j  is  the  penetrator  type)  centered  on  the  crossing  point  x  of  that  penetrator  (x  is  a  random 
variable).  Any  barrier  element  that  lies  within  the  detectability  width  of  a  given  penetrator 
(in  a  certain  realization  of  the  process)  is  said  to  be  detectable  by  that  penetrator  or 
vulnerab'e  to  that  penetrator,  or  we  say  that  that  penetrator  can  potentially  detect  such  a 
barrier  element.  If  penetrator  p  can  potentially  detect  barrier  element  {3,  p  actually  detects  J3 
with  probability  dj. 

C.  DERIVATION  OF  THE  DESIRED  QUANTITIES 

This  section  presents  the  probabilistic  arguments  used  to  derive  formulas  for  the 
three  basic  quantities  of  interest:  the  expected  number  of  barrier  elements  killed,  the 
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expected  number  of  shots  fired  at  a  given  barrier  element,  and  (where  possible)  the 
probability  distribution  of  the  number  of  shots  fired  at  a  given  barrier  element.  The 
algebraic  evaluation  of  the  final  formulas  is  presented  in  Appendix  B;  the  formulas  are 
presented  as  theorems  in  Section  E,  below. 

The  derivation  proceeds  in  two  parts.  Subsection  1  considers  one  specific  target 
(barrier  element  (3)  and  one  specific  shooter  (penetrator  p,  assumed  to  be  of  some  type  j). 
Subsection  2  integrates  the  results  of  Subsection  1  into  formulas  for  the  overall  quantities 
of  interest. 

1.  One  Penetrator,  One  Barrier  Element 

Consider  one  specific  barrier  element  p,  and  one  specific  penetrator  p.  We  assume 
that  p  is  of  type  j.  There  is  only  one  type  of  barrier  element,  and  we  are  tacitly  assuming 
that  the  results  here  do  not  depend  on  the  specific  barrier  element  considered.  This  is  true 
in  the  absence  of  edge  effects.  For  this  section  we  ignore  edge  effects;  some  of  the 
complications  they  can  cause  are  discussed  in  Section  D.  Unlike  Chapter  I,  we  do  not 
assume  that  P  is  detectable  by  p.  Instead  we  explicitly  consider  the  probability  of 
vulnerability. 

We  define  the  following  events: 

V— p  can  potentially  detect  P,  i.e.,  p's  crossing  point  is  within  wj/2  of  the  location 
of  p. 

Wm— The  total  number  of  barrier  elements  that  p  can  potentially  detect  is  exactly  m. 
(Defined  for  m=0,  1,...,B.) 

VWm--The  intersection  of  V  and  Wm,  i.e.,  p  can  potentially  detect  P  and  m-1 
other  barrier  elements. 

D--p  actually  detects  p. 

A^-p  detects  exactly  t  targets  other  than  target  p. 

DAi— The  intersection  of  D  and  Ai,  i.e.,  p  detects  target  P  and  exactly  i  other 
targets,  for  a  total  of  i+l.  (Defined  for  i=0,...,B-l.) 

Fh--p  fires  exactly  h  shots  at  p.  (Defined  for  h=0,...,Hj.) 

K-p  kills  p. 
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We  also  let  the  random  variable  K  denote  the  number  of  shots  p  fires  at  p,  so 

P(Fh)  =  P(K=h). 

The  quantities  we  derive  here  are  P(K),  E(K),  and  the  P(Fh)  for  h=0,...,H;.  These 
will  be  used  as  inputs  to  Subsection  2. 

By  Assumption  5),  target  (3  will  receive  no  more  than  Hj  shots  from  shooter  p  (in 
attack  protocols  (a)  and  (e),  Hj=l  for  all  j)  so  by  Assumption  7) 

P(K)  =  X  (l-fl-k.^P(Fh) 
h=0 

=  X  (1— [1— kj  l^1)  P(Fh  ). 
h=l 

Similarly 

C  VJ 

E(fT)  =  X  hp(ph>- 
h=l 


In  attack  protocols  (b),  (c),  (d),  and  (e),  the  number  of  shots  p  fires  at  P  depends  on  the 
total  number  of  targets  p  has  detected,  or  equivalently,  the  number  t  of  targets  other  than  P 
that  p  has  detected.  Assumption  5)  of  Chapter  II  is  identical  to  Assumption  5)  of  Chapter 
I,  except  that  penetrators  and  barrier  elements  have  reversed  shooter  and  target  roles.  Thus 
the  probabilities 

P(FhlDAi) 

for  h=0,...,Hj  and  the  appropriate  values  of  i  are  identical  to  those  used  in  Chapter  I,  for 
each  attack  protocol.  These  probabilities  are  presented  in  Lemma  1.5  of  Appendix  A.  (To 

use  these  in  context  here,  substitute  Hj  for  H.)  Since  D  and  {DAl}v>o  form  a  partition, 
and  a  shooter  fires  zero  shots  at  a  target  it  does  not  detect,  by  the  law  of  total  probability 

P(Fq)  =  P(D)  +  X  PCFoiDAJPfDAO 

l 
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and 


P(Fh)  =  £  P^lDAOPtDAO,  h=l,...,H. 

I 

To  deri”^  the  probabilities  P(DAi),  we  first  determine  the  number  of  targets  that  p  can 
potentially  detect  considering  the  spacing  of  the  targets  along  the  barrier  line.  Clearly, 
ignoring  edge  effects,  P(V)  =  wj/L  because  p's  crossing  point  is  uniformly  distributed  on 
[0,L]  (Assumptions  2)  and  4)).  We  note  that  (by  Assumption  4))  P(D)  =  djP(V)  = 
wjdj/L.  The  probabilities  of  the  Wm  (m=0,...,B)  can  be  computed  by  methods  similar  to 
those  used  in  Chapter  I,  but  V  and  Wm  are,  in  general,  not  independent.  The  following 
lemma,  which  is  of  key  importance  in  the  derivation,  is  proved  in  Appendix  B. 


Lemma  II.  1  Let  Ij  =  LwjB/lJ  and  fj  =  <WjB/L>.  Then,  ignoring  edge  effects, 

ijd-9 
P(W',IV>  ■  <9vT) 

dj+Dfj 


p(Wmlv)  =  0  for  all  other  values  of  m. 


(Parameter  L  is  restricted  to  be  >  0.  If  wj=0,  then  penetrator  p  cannot  potentially  detect  any 
target,  and  the  joint  probabilities  P(VWm)  are  zero  for  all  m.)  We  also  state  the 
probabilities  of  the  Wm. 

Lemma  II.2  Let  Ij  and  fj  be  defined  as  in  Lemma  II.l.  Then,  ignoring  edge 

effects, 

P(W,  )  =  1-fj 

P(Wi'+l>=fj 

p(Wm)  =  0  for  all  other  values  of  m. 


Lemma  II.2  is  proved  in  Appendix  B.  The  problem  of  edge  effects  is  discussed  in  Section 
D,  below.  The  joint  probabilities  P(VWj  )  and  P(VWj  +  j )  are  then 


P(VW,  )  =  PCWjIVjPCV)  = 


B 


n-6 


and 


(Ij+Dfj 

P(VWI+1)  =  P(W,.+1IV)  P(V)  = 

Suppose  that  Y\Ym  occurs.  Then  the  event  DAL  means  that  p  detects  P  and  exactly 
i  of  the  m-1  other  targets  vulnerable  to  p.  By  Assumption  4),  detections  of  different 
vulnerable  targets  are  independent  and  each  occurs  with  probability  dj.  Therefore  the 
number  of  other  targets  detected  follows  a  binomial  distribution  with  parameters  m-1  and 
dj,  and 


P(DAlIVWm)  =  ^ 


( 


\ 


m-1 

i 


j 


/i  j  jn-l-i 
dj  (1-Kip 


1=0,...,  m-1. 


Since  the  event  VWm  has  positive  probability  only  if  m=Ij  or  Ij  +1,  unconditionally, 
P(DA  )  =  P(DA  IVWT  )  P(VWT  1  +  P(DA  I  VWt  )  P(VWT  J_1  ). 

i  i  j  j  i  li+l  J  1 

Substituting  the  appropriate  derived  expressions  and  simplifying  yields 


Ui-fj) 
P(DA  )  =  — — ^  J— 


V' 

v  i 


i+l  I  -i-i 
^  ( 1— dj ) 


dj+1)fj 


B 


i  j 


d; 


l+l 


(1— dj) 


Is  -i 


1=0,...,^  -1 


and 


P(DAj  )  = 


(Ij+Dfj  lj+i 

-J - l  d- J 


If  lj=0,  then  the  first  formula  above  does  not  make  sense.  But  in  this  case,  at  most 
one  target  can  be  vulnerable  to  penetrator  p.  The  occurrence  of  D  implies  that  target  (5  is 
vulnerable  to  p,  thus  P  is  that  one  target  and  there  are  no  other  targets  that  p  can  potentially 

detect.  Thus  if  lj=0,  DA;  is  the  null  event  for  i  >  1,  and  DAo  is  simply  D,  so  P(DAo)  = 
wjdj/L. 

Note  also  that  if  Ij  >  0  and  dj=l,  P(DA;)  will  be  zero  for  all  i  except  Ij— 1  and  Ij. 
This  results  in  considerable  simplification  of  several  formulas. 

Overall,  for  h  >1 
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I 

J 

p(Fh)  =X  P(FhIDA  )P(DA  ). 
1=0  1  1 


Substituting  this  expression  into  the  expressions  for  P(K)  and  E(K)  we  obtain 

Hi  3, 

P(K)=1  X  (l-[l-kj]h)P(FhIDA)P(DA) 

u_i  J  V  V 

and 


•h=l 


i=0 


Hi  F 
J  J 

E(h)  =  II  hP(FhIDAt)  P(DAt). 

h=l  i=0 

As  in  Chapter  I,  these  expressions  are  most  easily  evaluated  by  interchanging  the  order  of 
summation.  The  quantities 

H, 

P(KIDA  )=£  (Hl-kjrt  P(FhIDA) 

1  h=l  1 

and 

E(friDA)  hPfFhlDAO 

1  h=l 

(the  lower  limits  of  summation  could  be  zero  as  well  as  one)  are  exactly  the  same  as  the 
corresponding  quantities  of  Chapter  I.1  Expressions  for  them  (for  each  different  attack 
protocol)  appear  in  Table  A-l  of  Appendix  A.  Then 

P(K)  =  X  P(KIDA  )  P(DA  ) 

£6  1  1 

and 

E(K)  =  5*  E(hlDA)P(DA), 

i  i 

i=0 

and  these  quantities  are  evaluated  simply  by  substituting  the  appropriate  previously  derived 
expressions  (noting  the  t  where  PfDA^)  =  0)  and  simplifying.  Results  for  each  protocol 

are  presented  in  Appendix  B. 


The  formulas  are  the  same  algebraically.  However,  parameters  wj,  kj,  and  dj  now  apply  to  penetrators, 
not  barrier  elements,  as  shooters,  and  will  in  general  have  different  numerical  values  from  the 
corresponding  Chapter  I  parameters,  so  the  resultant  values  will  differ.  Also,  Hi  should  be  substituted  for 
H. 
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2.  All  Penetrators  and  Barrier  Elements 

To  use  the  results  of  the  previous  section  to  find  the  overall  quantities  of 
interest,  we  introduce  new  symbols  for  P(K),  E(K),  and  P(K=h),  as  follows.  For  j=l,...,n 
define 

qj--probability  a  given  barrier  element  is  killed  by  a  given  penetrator  of  type  j. 
(Ignoring  edge  effects,  this  probability  is  the  same  for  all  barrier  elements.) 

ej— expected  number  of  shots  a  given  penetrator  of  type  j  fires  at  a  given  barrier 
element. 

4>hj-probability  a  given  penetrator  of  type  j  fires  exactly  h  shots  at  a  given  barrier 
element.  (Defined  for  h=0,...,Hj.) 

Also  define: 

A  B-Number  of  barrier  elements  killed 
E(AB)— Expected  number  of  barrier  elements  killed 

S-Event  that  a  given  barrier  element  survives  every  penetrator 
y-Number  of  shots  fired  at  a  given  barrier  element,  totaled  over  all 
penetrators. 


We  seek  expressions  for  E(AB),  E(y),  and  P(y=y)  fo»  appropriate  values  of  y.  We  present 
the  main  issues  of  derivation  for  each  quantity  in  turn. 

By  elementary  arguments  [12], 


E(AB)  =  BP(S)  =  B(1-P(S)). 

The  event  S  is  the  intersection  of  the  events  that  the  given  barrier  element  survives  a  given 
penetrator,  taken  over  all  penetrators.  By  Assumptions  3)  and  6),  these  individual  survival 
events  are  independent.  Since  there  are  Rj  penetrators  of  type  j,  for  j=l,..„  n, 

P(S)  =  n  (l-q.)  J 
j=l  J  • 


Thus 


E(AB) 


B 


n  R. 

l-nO-^d 
j=l  J 
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Similarly,  y  is  the  total  over  all  penetrators  of  the  random  variables  Kp,  the  number 
of  shots  fired  at  the  given  barrier  element  by  penetrator  p  (p=l,...,R).  If  p  is  of  type  j, 

E(hp)  =  e^,  as  defined  above.  Therefore 

E(y)  =  tv,  • 

J=1 

In  fact,  by  Assumptions  3)  and  6),  the  random  variables  hp,  for  different  p  are 
mutually  independent.  In  attack  protocols  (a),  (d),  and  (e),  each  hp  can  assume  only  the 
values  zero  and  one,  and  if  p  is  of  type  j 

P(hp=l)  =  0lj  , 

and  the  <bij  can  be  different  for  different  j.  If  there  is  only  one  type  of  penetrator,  y=^Thp 

follows  a  binomial  distribution  with  parameters  R  (i.e.,  Ri)  and  <|>n.  If  n  >  1,  we  can  use 
the  Poisson  approximation  to  the  binomial  and  say  that  y  is  approximately  distributed 
Poisson  with  mean 


n 

H  -  X  RjV 

pl 

The  above  discussion  applies  only  if  attack  protocol  (a),  (d),  or  (e)  is  used.  If  attack 
protocol  (b)  or  (c)  is  used,  each  hp  can  assume  a  range  of  values,  and  even  though  the  hp 
are  independent  for  different  p,  we  have  not  obtained  the  resulting  distribution  for  y.  One 
simplifying  case  occurs  if  wj  <  L/B  for  all  penetrator  types  j.  In  this  case,  each  penetrator 
can  potentially  detect  no  more  than  one  barrier  element,  and  thus  will  fire  all  of  its  Hj  shots 
at  any  barrier  element  it  detects.  Thus  if  p  is  of  type  j,  Kp  can  assume  only  the  values  0  or 

Hj.  If  in  addition  there  is  only  one  type  of  penetrator,  then  y,  the  number  of  shots  fired  at  a 
given  barrier  element,  equals  Hi  multiplied  by  the  number  of  penetrators  that  detect  that 
barrier  element.  If  I  j  =0,  this  latter  quantity  is  a  binomial  random  variable  with  parameters  R 
(i.e.,  Rj)  and  <J>Hj  i-  Therefore,  if  n=l,  1=0,  and  attack  protocol  (b)  or  (c)  is  used. 


All  of  the  above  formulas  were  expressed  in  terms  of  the  qj,  ej,  and  <|>hj.  The 
algebraic  expressions  derived  in  Subsection  1  for  P(K),  E(fT),  and  P(h=h)  can  be 


11-10 


substituted  in  and  overall  algebraic  formulas  computed.  The  final  results  are  presented  in 
Section  E;  appropriate  details  appear  in  Appendix  B. 


D.  EDGE  EFFECTS 

Edge  effects  will  occur  if  some  barrier  elements  are  located  less  than  wj/2  away 
from  the  edge  of  the  barrier,  for  some  penetrator  type  j— i.e.,  if  wj  >  L/B.  This  condition  is 
similar  to  the  condition  in  Chapter  1,  Section  F,  but  wj  is  now  a  penetrator' s  detectability 
width.  This  has  two  consequences,  which  work  in  opposing  directions.  Consider  a 
specific  barrier  element  (3  located  less  than  wj/2  away  from  the  barrier  edge,  for  some  j,  and 
consider  a  penetrator  p  of  type  j.  The  probability  P(V)  that  p  can  potentially  detect  P  is 
lower  than  wj/L,  which  was  used  in  the  computations  in  Section  C.  However,  the  total 
number  of  barrier  elements  p  can  potentially  detect  assumes  values  lower  than  Ij  or  Ij  +  1  if 
p's  "detectability  width"  extends  beyond  the  barrier  edge.  If  an  attack  protocol  is  used  in 
which  an  attacker  has  a  limited  number  of  shots,  if  p  does  detect  p  there  are  fewer  other 
detectees  to  divert  p's  shots,  and  the  probability  p  is  killed  increases. 


If  edge  effects  are  considered,  different  penetrators  of  the  same  type  can  still  be 
considered  as  identical  (and  mutually  independent),  but  different  barrier  elements  cannot. 
The  probability  qj  that  a  specific  penetrator  of  type  j  kills  a  specific  barrier  element  must  be 
replaced  by  qjb,  the  probability  that  the  penetrator  kills  barrier  element  b,  for  b=l,...,B. 


Then 


Similar  arguments  apply  to  the  expectation  and  probability  distribution  of  y  . 


If  wj  <  L/B,  then  edge  effects  do  not  occur  for  penetrators  of  type  j,  and  the 
probabilities  P(VWm)  are  as  given  in  Section  C.  The  reason  for  this  is  discussed  in 
Appendix  B.  Since  the  overall  quantities  of  interest  involve  sums  or  products  on  j,  wj  must 
be  less  than  or  equal  to  L/B  for  all  j  if  edge  effects  are  not  to  occur.  (Unlike  Chapter  I,  this 
condition  is  not  needed  to  guarantee  the  independence  of  any  set  of  events.)  The  phrase 
"ignoring  edge  effects"  in  the  statements  of  the  theorems  in  the  next  section  merely  means 
that  if  wj  <  L/B  for  all  j,  the  formulas  are  exact;  if  not,  the  formulas  are  valid  except  for  edge 
effects. 

We  note  finally  that  although  the  combat  processes  examined  in  this  paper 
considered  barrier  elements  spaced  along  a  line,  one  could  also  consider  a  circular  barrier  or 
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screen,  with  barrier  elements  on  the  circumference  (which  is  of  length  L).  If  detectability 
areas  are  considered  as  portions  of  the  circumference,  the  results  in  both  Chapters  I  and  II 
hold  without  edge  effects  regardless  of  the  relations  between  the  wj  and  L/B. 

E.  THE  RESULTS 

This  section  presents  the  formulas  for  the  three  main  quantities  of  interest:  E(/£*B), 
the  expected  number  of  barrier  elements  killed;  E(y),  the  expected  number  of  shots  fired  at  a 
given  barrier  element;  and  the  probabilities  P(y=y)  that  exactly  y  shots  are  fired  at  a  given 
barrier  element,  for  appropriate  values  y.  (The  P(y=y)  have  been  derived  in  general  only  for 
attack  protocols  (a),  (d),  and  (e).)  The  methods  of  derivation  of  the  formulas  have  been 
presented  in  Section  C;  proof  details  appear  in  Appendix  B.  Results  for  which  we  have  not 
found  a  closed  form  are  expressed  with  cumulative  binomials,  when  possible. 

The  meanings  of  the  parameters  B,  R,  n,  L,  wj,  dj,  Hj,  and  kj  are  given  in  Section 
A,  above.  We  let  Ij  and  fi  denote  the  integer  and  fractional  parts,  respectively,  of  wjB/L. 
Note  that  Ij  =  0  if,  and  only  if,  wj  <  L/B.  It  can  easily  be  verified,  however,  that  the 
formulas  given  below  for  the  case  where  Ij  =  0  also  hold  if  wj  =  L/B,  i.e.,  Ij  =  1  and  fj  =  0. 
For  those  penetrator  types  j  where  Ij  >  1,  we  define  the  parameters 

+  F(Ij  +l)fJdj  i'/k  for  i=G,...,Ij-l, 

and 

V  B<Ii+1)fid/,+  1. 

(Note  timi  0j t  is  the  probability  P(DAt)  used  in  Section  C.)  As  in  Chapter  I,  Section  G,  we 
use  the  following  notation  for  binomial  probabilities: 

b(m  ;  n,p)  =  f "  )pmG-p)n"m  m=0,...,n 

and 
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B(m;n,p) 


r 


0, 


b(i ;  n,p). 


m  <-l 
m=0,...,n 


1, 


m  >  n. 


The  theorems  are  indexed  by  attack  protocol,  reflecting  dependence  of  the  results  on 
the  attack  protocol.  The  expected  number  of  penetrators  killed  is  given  by  the  following 
theorems. 


Theorems  Tl.la.  b.  e.  d.  e.  Ignoring  edge  effects,  the  expected  number  of 
penetrators  killed  is  given  by 

E(Afe)  =  B  [1  -  FI  (l-qj)Rjl  , 
j=1 

where  the  qj  are  defined  separately  for  each  attack  protocol  (a)  through  (e),  as  follows. 

-  Attack  protocol  (a) 

q;  =  kjWjdj/L 

-  Attack  protocol  (b) 

If  Hj  =  0,  q.  =  0. 

If  Hj  >  1  and  I.  =  0, 


If  Hj  >  1  and  Ij  >  1 


-  Attack  protocol  (cl 
If  Hj  =  0,  qj  =  0. 

If  ^  >  1  and  Ij  =  0, 


as  in  protocol  (b). 


Hi  w-dj 

qj  =  [1-  (1-k  )  . 

J  1— / 


W;  d; 


k.  H 


‘j  L  l+l 


qj  =  [1  -  d-kj)  J  ] 


H:  Wjdj 
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If  Hj  >  1  and  L  >  1 


!j 


H. 


qj  =  in  -  o-kj)L 


l+l 


(1  - 


1=0 


<lT7>kj>)  ei 


-  Attack  protocol  (dl 


IfHj  =  0,qj  =  0. 

If  1  <  Hj  <  Ij 

ij-  Tho-^i  B<Hj-|;'j-1-dj) 

+  ^<Ii+1>fjdj®(Hi-l;lj.  dj) 

kj 

+  [1  -  B(Hj;  \f  dp  ] 

k, 

+  ^  H.f.  [1  -  B(Hj;l.+l,  dj)]  . 
if  ij  =  0  or  Hj  *  !j  + *>  qj  =  kJwjdj/L  • 

-  Attack  protocol  (el 

k-  j. 

qj=_i_n  . 

If  Ij  =  0,  this  reduces  to 

Ij  *  ^Wjdj/L. 

The  above  results  for  attack  protocols  (b)  and  (c)  do  not,  in  general,  have  closed 
forms.  If  the  detection  probability  dj  equals  1,  however,  most  of  the  0ji  are  zero  and  the 
formulas  for  qj  are  easily  evaluated.  In  the  implementation  of  these  protocols  in  the 
NAVMOD  model  [3],  parameters  wj  and  dj  are  input  but  the  code  evaluates  the  expression 
for  qj  using  wjdj  and  1  instead  of  wj  and  dj. 

It  can  be  verified  that  if  Hj  is  equal  to  1 ,  the  formulas  for  qj  in  attack  protocols  (b), 
(c),  and  (d)  all  reduce  to  the  corresponding  formulas  in  protocol  (e). 
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The  expected  number  of  shots  fired  at  a  specific  barrier  element  is  given  by 

Theorems  II.2a.  h.  c.  d.  e:  Ignoring  edge  effects,  the  expectation  E(y)  of 
the  number  of  shots  y  fired  at  a  particular  barrier  element  (5  by  the  penetrators  is 

n 

E(y)  =  jRjej  , 
j=l 

where  the  ej  are  defined  separately  for  each  attack  protocol,  as  follows. 

-  Attack  protocol  (a) 

ei  ■  wid,/L  • 

-  Attack  protocol  (b) 

ej-£n  -a-d/a-WJ  • 

-  Attack  protocol  (cl 

'i-T11  -(1VJ(1-W1  • 

-  Attack  protocol  (dl 

ej  -  ij*j  • 

where  qj  is  as  defined  in  Theorem  II. Id,  for  the  appropriate  combination  of  Hj  and  Ij 
values. 

-  Attack  protocol  (el 

ei  -  F  t‘  -  <‘V’  <‘-fj A  >  1 

In  accordance  with  the  arguments  of  Section  C.2,  the  probability  distribution 
of  y ,  the  total  number  of  shots  fired  at  a  given  barrier  element,  is  given  for  the  following 
cases. 

-  There  is  only  one  penetrator  type  (n=l)  and  attack  protocol  (a),  (d),  or  (e)  is  used. 

-  There  is  only  one  penetrator  type,  attack  protocol  (b)  or  (c)  is  used;  and  wj  <  L/B. 

-  There  are  multiple  penetrator  types  and  attack  protocol  (a),  (d),  or  (e)  is  used. 

Theorems  II.3a.  d.  e:  Let  n=l.  Then,  ignoring  edge  effects,  y,  the  number  of 
shots  fired  at  a  particular  barrier  element,  is  binomially  distributed  with  parameters  R  and 
<J>1  j,  where  <J>i  \  is  defined  separately  for  attack  protocols  (a),  (d),  and  (e),  as  follows. 
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-  Attack  protocol  (a) 


4>n  =  widi/L 

-  Attack  protocol  (di 

If  Hi  =  0,  <t>n  =  0  . 

If  Hj  >  1  and  Ij  =  0 

0ii  =  Wjdj/L 

If  Hj  >  1  and  It  >  1 

^  i!  (1-fj )  dj  B(Hj  -1;  Ij-1,  dj ) 

+  i-  (I,  +1)  fj  dj  I! ,  dj ) 

+  1  Hj  (1-fj)  [1  -  ®(Hj ;  Ilt  dj )  ] 

+  ^  H,fj  [1  -B(H;  Ii+1,  dj)]. 

(This  is  identical  to  e  ]  and  qj/kj.)  It  can  be  -fied  that  if  Hj  >  Ij  +  1,  the  above 
formula  reduces  to 

$11  =  Wjdj/L  . 

-  Attack  prQtocpl  (gJ 

$n  =  B  [1  — (1— di  Z1  (l_fl  di  )  1  • 

The  subscripts  on  (J)^  are  not  really  necessary;  we  have  kept  them  because  <t>n  is  the  same 

quantity  indicated  by  that  symbol  in  Section  II. C,  above:  the  probability  a  given  penetrator 
fires  exactly  one  shot  at  the  particular  barrier  element  considered. 
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and  P(y=y)  =  0  for  all  other  values  of  y.  (I.e.,  y/Hj  is  binomially  distributed  with 
parameters  R  and  widi/L.) 


Suppose  that  attack  protocol  (a),  (d),  or  (e)  is  used. 


Then,  possibly  subject  to  edge  effects,  y  is  approximately  Poisson  distributed  with  mean 


where  the  ej  are  as  given  in  Theorems  II.2a,  d,  and  e  above. 


(Since  in  attack  protocols  (a),  (d),  and  (e),  a  shooter  fires  no  more  than  one  shot  at  any 
particular  target,  the  quantity  ej,  which  is  the  expected  number  of  shots  fired  at  a  given  target 
by  a  given  shooter,  is  the  same  as  the  probability  the  shooter  fires  one  shot  at  that  target.) 
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PROOFS  OF  RESULTS  IN  CHAPTER  I 


This  appendix  is  in  two  parts.  The  first  part  proves  seven  lemmas  that  are  used  in 
deriving  the  overall  formulas.  The  first  four  were  stated  in  Chapter  I,  the  others  provide 
useful  algebraic  results.  The  lemmas  involve  the  following  issues: 

Lemma  1. 1— Independence  condition  for  attack  protocol  (a). 

Lemma  I.2-Independence  condition  if  wj  <  L/B. 

Lemma  I.3--Formula  for  Nj(x). 

Lemma  1.4— Formula  for  P(y=y),  in  terms  of  P(Fy)  when  wj  <  L/B. 

Lemma  1.5— Derivation  of  probabilities  P(FhlDAt). 

Lemma  1.6-Formulas  for  P(K|DAt)  and  E(KlDAv)  . 

Lemma  1.7— Formulas  for  sums  and  weighted  sums  of  (Ox  and  Yji- 

The  second  part  of  the  appendix  presents  the  steps  in  proving  the  theorems  of 
Chapter  I,  Section  G.  The  main  conditioning  arguments  have  been  presented  in  Chapter  I. 
Section  E.  Here,  these  will  be  briefly  restated,  with  additional  arguments  as  necessary  to 
complete  the  proofs. 

Unless  specifically  redefined,  the  notation  used  in  this  appendix  matches  that  of 
Chapter  I.  The  parameters  Rj,  R,  B,  n,  L,  wj,  dj,  kj,  and  H,  and  the  attack  protocols  (a) 
through  (e)  are  as  defined  in  Chapter  I,  Section  A.  We  let  Ij  and  fj  denote  the  integer  and 
fractional  parts,  respectively,  of  WjB/L.  The  events  and  random  variables  Fh,  D,  DAl(  K, 
and  Ti  concerning  one  penetrator  and  one  barrier  element  are  as  defined  in  Chapter  I, 
Section  E.l.  The  events  and  random  variables  S,  Sb,  Cx,  x,  Nj(x),  >’b»  and  y  are  as 

defined  in  Chapter  I,  Section  E.2. 
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PROOFS  OF  RESULTS  IN  CHAPTER  I 


This  appendix  is  in  two  parts.  The  first  part  proves  seven  lemmas  that  are  used  in 
deriving  the  overall  formulas.  The  first  four  were  stated  in  Chapter  I,  the  others  provide 
useful  algebraic  results.  The  lemmas  involve  the  following  issues: 

Lemma  I.i— Independence  condition  for  attack  protocol  (a). 

Lemma  I.2-Independence  condition  if  wj  <  L/B. 

Lemma  I.3--Formula  for  Nj(x). 

Lemma  1.4— Formula  for  P(y=y),  in  terms  of  P(Fy)  when  wj  <  L/B. 

Lemma  1.5— Derivation  of  probabilities  P(FhlDAi). 
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Lemma  1.7— Formulas  for  sums  and  weighted  sums  of  g>i  and  \j/ji. 

The  second  pan  of  the  appendix  presents  the  steps  in  proving  the  theorems  of 
Chapter  I,  Section  G.  The  main  conditioning  arguments  have  been  presented  in  Chapter  I, 
Section  E.  Here,  these  will  be  briefly  restated,  with  additional  arguments  as  necessary  to 
complete  the  proofs. 

Unless  specifically  redefined,  the  notation  used  in  this  appendix  matches  that  of 
Chapter  I.  The  parameters  Rj,  R,  B,  n,  L,  wj,  dj,  kj,  and  H,  and  the  attack  protocols  (a) 
through  (e)  are  as  defined  in  Chapter  I,  Section  A.  We  let  Ij  and  fj  denote  the  integer  and 
fractional  parts,  respectively,  of  wjB/L.  The  events  and  random  variables  Fh,  D,  DAV,  K, 
and  H  concerning  one  penetrator  and  one  barrier  element  are  as  defined  in  Chapter  I, 
Section  E.l,  The  events  and  random  variables  S,  Sb,  Cx,  x,  Nj(x),  yb*  anc*  y  are  as 

defined  in  Chapter  I,  Section  E.2. 


A.  THE  LEMMAS 


Lemma  1.1  Consider  one  penetrator  p,  of  type  j.  Let  the  events  Sb,  S,  and  Cx  be 
as  defined  in  Chapter  I,  Section  E.2.  Suppose  that  attack  protocol  (a)  is  used.  Then,  for 
every  x  €  [0,  L] 

p(sicx)  =  n  p(sb'cx) . 

b=l 

Proof:  Suppose  that  Nj(x)  barrier  elements  can  potentially  detect  the  penetrator  if 
fbe  penetrator  crosses  at  point  x.  For  now,  let  N  denote  Nj(x).  By  Assumption  4),  the 
number  m  of  barrier  elements  that  actually  detect  the  penetrator  is  binomially  distributed 
with  parameters  N  and  dj.  Since  attack  protocol  (a)  is  used,  m  is  also  the  number  of  shots 
fired  at  penetrator  p.  By  Assumption  7),  the  probabilitv  p  survives,  given  that  m  shots  are 
fired  at  it,  is  (l-kj)m,  regardless  of  p's  crossing  point.  Then 


Sill  A  ^  111 

( 1  — k j )  m  dj  (1— dj ) 

m-0  v  J 


Using  the  formula  for  the  z-transform  of  the  binomial1  this  becomes 


PCS  ICX )  =  [l-dj+dj(l-kj)]N  =  (l-djkj)^00. 

Now  consider  barrier  element  b.  Whether  b  can  detect  p  depends  on  p's  crossing 
point  x.  Given  x,  if  b  cannot  potentially  detect  p,  p  will  certainly  survive  b  and 

P(SbICx)  =  1. 

If  b  can  potentially  detect  p,  then  b  will  actually  detect  p  with  probability  dj.  Since  protocol 
(a)  is  used,  b  will  then  fire  one  shot  at  p,  which  will  be  lethal  with  probability  kj.  Thus  if  b 
can  potentially  detect  a  penetrator  crossing  at  x 

P(SbICx)  =  1  -  djkj  . 

Since  given  Cx,  Nj(x)  penetrators  can  potentially  detect  p, 


n 

m  n 

2JL  b(m;n,p)  =  (1-p+pz) 

m=0 

where  the  notation  of  Chapter  I,  Section  G  for  binomial  probabilities  has  been  used. 
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B  N  (x) 

n  p(sbicx)  =  y  jv  ; 

b=i  J  J 

which  equals  the  previously  derived  formula  for  P(SICX). 

Lemma  1.2  If  penetrator  p  is  of  type  j  and  wj  <  L/B,  then  the  events  Sb  that  p 
survives  barrier  elements  b,  for  b=l,...,B,  are  mutually  independent  given  the  crossing 
point  x  of  p. 

Proof:  Let  the  events  S  and  Cx  be  as  defined  in  Chapter  I,  Section  E.2.  By 
Lemma  1.3,  below,  the  number  of  barrier  elements  that  can  potentially  detect  p  is  Ij  with 
probability  1  -  fj  and  Ij  +  1  with  probability  fj.  If  wj  <  L/B,  either  Ij  =  0  or  Ij  =  1  and  fj  = 
0,  so  that  with  probability  one  p  is  vulnerable  to  at  most  one  barrier  element.  If  x  is  such 
that  p  is  vulnerable  to  no  barrier  element,  then  p  survives  each  barrier  element,  so 

P(SICX)  =  1  =  n  P(SbICx). 

If  x  is  such  that  p  is  vulnerable  to  one  barrier  element  P,  then  p  survives  overall  if  and  only 
if  it  survives  P- so  P(SplCx)  =  P(SICX).  For  b  *  P,  p  is  not  vulnerable  to  barrier  element 
b,  thus  P(St>ICx)  =  1.  Thus 

B 

P(SICX)  =  n  P(SbICx) 

b=i  D  * 

as  stated. 


Lemma  1.3  Let  the  random  variable  x  be  the  crossing  point  of  some  specific 
penetrator  p  of  type  j.  Let  Nj(x)  be  the  number  of  barrier  elements  that  can  potentially 
detect  that  penetrator.  Then,  ignoring  edge  effects, 


N.(x) 


I. 


J 

I.  +  1 
J 


w.p  1-f. 

r  J 
w.p.  f 


w.B 

and  f.  =  -j - I.  .  Alternatively  stated,  the  set  of  x  e  [0,L]  where 

J  L  J 

Nj(x)  =  Ij  has  length  (l-fj)L,  and  the  set  of  x  6  [0,L]  where  Nj(x)  =  Ij  +  1  has  length  fjL. 

Proof:  Let  the  barrier  line  (of  length  L,  by  Assumption  1))  be  labeled  with 
coordinates  0  through  L.  The  number  of  barrier  elements,  B,  is  assumed  to  be  an  integer  > 
l .  By  the  positioning  described  in  Assumption  1),  the  coordinate  of  barrier  element  b  is 


where  I.  = 
J 


w.B 

J 
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L  1. 
xb  =  g  (b-y)  • 

By  Assumption  2),  x  is  distributed  uniformly  on  [0,L].  For  any  specific  crossing  point 
x  e  [0,L]  of  penetrator  p,  by  Assumption  4),  barrier  element  b  will  detect  penetrator  p  if 
and  only  if 

lxb-xl  <Wj/2  . 


Substituting  ^  (b  _-I)  for  xb,  denoting  wj/L  by  8,  and  performing  the  algebra  yields  the 
equivalent  condition 


y  +  B(r 


)  <  b  <  j  +  B  +  B9 


It  is  clear  that  Nj(x)  is  the  cardinality  of  the  set  of  b  that  satisfy  the  above  condition.  We 
ignore  edge  effects  by  letting  b  assume  any  integer  values,  possibly  nonpositive  or 
exceeding  B.  (That  is,  p  crosses  in  the  middle  of  an  infinitely  long  barrier.)  Denote  the 
left-hand  side  of  the  above  expression  by  bp(x),  i.e.,  for  each  x  e  [0,L]  let 


b0(x)  =  y+  B(i-|)  . 

Denote  (l-B0)/2  by  ap.  Then,  since  x  ~  U(0,L],  bp(x)  is  a  random  variable  uniformly 
distributed  on  [ap,  ap+B],  and  is  integer  with  probability  zero.  The  desired  quantity  Nj(x) 
is  the  number  of  integer  points  in  the  interval 


[b0  (x)  ,  b0(x)  +  B0] . 

To  determine  this  number,  let  us  simplify  the  notation.  Let  bp(x)  be  denoted  by  a.  For  any 
noninteger  a  and  any  nonnegative  y,  the  integer  points  in  the  interval  [a,  a+y]  are  LaJ  +1, 
LaJ  +  2,...,  up  to  LaJ  +  k,  where  k  is  the  greatest  integer  such  that 

LaJ  +  k  <  a  +  y  . 


Rearranging,  we  see  that  k,  which  is  the  number  of  integer  points  in  the  interval,  equals 

L  (a+y  -  LaJ  )  J  , 
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which  equals 


LyJ  +  L<a>  +  <y>  J  . 


Since  fractional  parts  <a>  and  <y>  are  both  in  the  interval  [0,1), 

,  ,  / 0  if  <a>  <  l-<y> 

L  <a>  +  <y>j  =  if  <a>  >  l-<y>  . 

Now  let  a  be  a  random  variable  distributed  uniformly  on  [ao,  ao  +  B].  Remember  that  B  is 
the  number  of  barrier  elements,  a  positive  integer;  y  is  a  nonnegative  constant.  Define  the 
random  variable 

k  =  LyJ  +  L  <a>  +  <y>  J  . 

The  interval  [ao,  ao,  +  B]  is  composed  of  B  unit  intervals.  As  a  ranges  over  the  interval 
[ao,  ao  +  1),  <a>  assumes  every  value  between  0  and  1  (except  for  1  itself).  The  same  is 
true  as  a  ranges  over 

[a0+ 1,  a0 +2),  [a0+2,  a0 +3),..,[aQ+ (B-l),  a0+B)  . 

Therefore  <a>  is  uniformly  distributed  on  [0,1]  and  is  thus  less  than  1  -  <y>  with 
probability  1  -  <y>.  Then  the  distribution  of  1c  is 

/  LyJ  w.p.  1  -  <y> 

£  =  < 

^  LyJ+  1  w.p.  <y>  . 

But  as  it  has  been  defined,  lc  is  in  fact  Nj(x).  Let  y  =  B0  =  wjB/L.  Substituting  these  in 
the  above  expression  yields  the  statement  of  the  lemma. 

We  note  that  if  wj  <  L/B,  then  for  every  x  e  [0,L]  the  condition 

I  +  B  (L  "  -  b  -  2  +  B(L  "  +  80 

for  barrier  elements  b  that  can  potentially  detect  a  penetrator  crossing  at  x  strictly  includes  at 
most  one  integer  value  of  b,  and  this  value  is  not  negative  or  strictly  greater  than  B.  This  is 
true  because  wj  <  L/B  if  and  only  if  B6  <  1 .  The  interval  in  which  b  can  lie,  therefore, 
has  width  <1.  Its  left  endpoint  is  greater  than  or  equal  to  Bx/L,  which  is  never  negative 
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for  x  €  [0,L],  Similarly,  the  right  endpoint  is  <  Bx/L  +  1,  so  the  interval  includes  the 
integer  point  B  +  1  only  if  x  =  L,  which  happens  with  probability  zero.  Thus,  if  wj  < 
L/B,  the  edge  effect  of  "detection"  of  a  penetrator  by  mythical  barrier  elements  will  not 
occur. 


Lemma  1.4  Suppose  wj  <  L/B  for  all  i  =  l,....,n.  Then  for  all  attack  protocols 

-  Ifwj  <  L/B 

P(y=0)  =  1  -  fj  +  fjPCFo) 

P(y=y)  =  fjP(Fy) ,  y=l,...,H, 

-  Ifw.  =  L/B 

J 

P(y=y)  =  P(Fy),  y=0,...,H, 


where  in  attack  protocols  (a)  and  (e)  H  is  taken  to  be  1. 

Proof;  Let  yj,  be  the  number  of  shots  fired  at  the  given  penetrator  p  (of  type  j)  by 
barrier  element  b,  so  y,  the  total  number  of  shots  fired  at  p,  is 


y  = 


B 


Fy  is  the  event  that  a  given  barrier  element  p  fires  exactly  y  shots  at  a  given  penetrator  p  (of 
type  j).  The  probabilities  P(Fy)  are  calculated  assuming  that  p  can  potentially  detect  p.  Let 
N  denote  the  number  of  barrier  elements  that  can  potentially  detect  p.  By  Lemma  1.3, 


w.p.  1-fj 

w.p.  f;  , 

and  if  wj  <  L/B,  edge  effects  are  not  a  problem.  A  barrier  element  that  cannot  potentially 
detect  p  fires  zero  shots  at  p.  Thus  P(y=0  I  N=0)  =  1  and  for  y=l,2,..„  P(y=y  I  N=0)  =  0. 
Also  P(y=y  I  N=l)  is  simply  P(Fy),  for  y=0,l,2,...,  because  if  N  =  1,  the  number  of 
shots  p  receives  overall  is  simply  the  number  of  shots  p  receives  from  the  one  barrier 
element  that  can  potentially  detect  p. 
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Always, 


P(y=y)  =  P(y=y  iN=lj  Xi-fj )  +  P(y=y  iN=ij  +D  fj  • 

If  wj  <  L/B,  Ij  =  0,  and  if  wj  =  L/B,  Ij  =  1  and  fj  =  0.  Making  the  appropriate 
substitutions  for  the  indicated  conditional  probabilities  in  the  above  equation  yields  the 
statements  in  the  lemma.  (Recall  from  Chapter  I,  Section  E  that  P(Fy)  =  0  for  y  >  H.) 

Lemma  1.5  Consider  one  attacker,  barrier  element  P,  and  one  target,  penetrator 
p.  Assume  that  {3  has  detected  p  and  exactly  i  other  targets.  Then  the  probability 
P(FhlDAt)  that  p  fires  exactly  h  shots  at  p  is  given  by  the  following  formulas,  for  each 
attack  protocol  in  turn  (i  can  be  any  nonnegative  integer). 


-  Attack  protocol  (a) 


PfFjIDAO  =  1 

P(Fh  I  DAt )  =  0  for  all  other  h. 


-  Attack  protocol  (b) 

P(FhIDAt)  = 


vS 


(-L-)h  (1_  -i-)H_h 


h=0 . ,H. 


i+l 


l+l 


H 


-  Attack  protocol  (c)  Let  J  denote  the  integer  part  and  Tj  the  fractional  part  of  •  Then 


PfFjIDAO  =  1  -  T\ 
P(FJ+1IDAV)  =  Ti 
P(FhIDA0  =  0  for  all  other  h. 

-  Attack  protocol  (d). 

P(F0iDA^)  =  1  -  min(l,  -—■) 

P(F,  IDA  )  =  min (1, -Hr-) 

1  i  v+1 

P(Fh IDA  )  =  0  h>2. 
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-  Attack  protocol  (e) 


P(F0 1  DA  )  =  1 


P(F,  i  DA  )  =  777- 


P(Fh  I  DA^)  =  0,  h  >2. 

Proof:  In  all  cases,  the  results  follow  from  the  statement  of  the  attack  protocols  in 
Assumption  5)  of  the  problem  statement  in  Chapter  I,  Section  A,  where  the  total  number  of 
targets  p  detects  is  1  +  1.  The  formulas  for  attack  protocols  (a)  and  (e)  are  self- 
explanatory.  The  formula  for  (b)  is  a  precise  statement  of  the  attack  protocol.  In  protocol 
(c),  each  target  will  receive  J  or  J  +  1  shots,  and  J  shots  at  each  target  uses  up  J(i+1) 
shots,  so  there  are  H  -  J(t+1)  shots  left  over,  to  be  distributed  among  the  1  +  1  targets 
such  that  no  target  gets  more  than  one  of  the  additional  shots.  The  probability  that  a  target 
does  get  one  of  the  additional  shots  can  be  seen  to  be 


-  J(i+D 

1  +  1 


-  J  =  T1  • 


For  attack  protocol  (d),  if  H  >  1  +  1,  then  each  target  will  receive  one  shot,  i.e., 
P(Fjl  DAO  =  1.  If  H  <  i+  1,  there  are  H  shots  to  be  apportioned  among  1  +  1  targets 
such  that  no  target  gets  more  than  one  shot.  The  probability  that  P  gets  one  of  the  shots  is 
simply  the  number  of  subsets  of  the  1  +  1  targets  of  size  H  that  contain  P,  divided  by  the 
total  number  of  subsets  of  size  H  of  the  1  +  1  targets,  i.e.. 


>)  =  [h-i)/(1h 


P(F,  IDA 


which  simplifies  to  — p  .  Combining  the  two  results,  we  obtain,  as  stated 


P(F  ,  I  DA  )  =  min  (1,  -2-)  . 

1  i+l 

Let  the  probabilities  P(Fh  I  DA^)  be  as  given  in  Lemma  1.5.  Then  the 


quantities 
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and 


P(KIDAt)  =  £  (l-[l-kj]h)P(FhIDAi) 

h=0 

H 

E(HlDAi)  =£hP(FhDAO 
h=0 

are  as  given  in  Table  A-l.  (Note  that  i  can  bu  any  non-negative  integer.  In  the 
homogeneous  case,  substitute  k  for  kj.) 

Proof:  The  only  issue  is  the  algebra,  which  is  straightforward  except  for  the  evaluation  of 
P(KIDA^)  in  attack  protocol  (b).  There  we  use  the  formula  for  the  z-transform  of  the 
binomial 

Y  zmb  (m  ;  n,  p)  =  (l-p+pz)n  . 
m=l) 

In  protocol  (c),  note  that  +  <i^T>  =  T+P  simP^^es  the  formula  for 

E(FTIDA^)  The  formulas  for  E(FT  IDA^)  are  identical  for  attack  protocols  (b)  and  (c).  This 

is  reasonable  given  the  nature  of  the  protocols.  In  (b),  an  attacker  scatters  its  H  shots 
randomly  over  the  t  +  1  targets  it  has  detected;  in  (c)  it  distributes  the  shots  as  evenly  as 

possible,  but  in  either  case  there  is  an  average  of  — shots  per  detected  target.  Since  the 

E(Ti  I  DA  )  are  the  same  for  all  t,  the  overall  E(h)  and  E(y)  are  the  same  for  protocols  (b) 
and  (c). 
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Table  A-l.  FORMULAS  FOR  PCKlDA^  AND  E(h  |DAt) 


Attack  Protocol 

PCKIDA^ 

E(hlDAt) 

(a) 

kj 

1 

(b) 

1 

1 

H 

i+l 

(c) 

H 

[i+lj 

l-(l-k.)  (i_<-H_>kj) 

H 

i+l 

(d) 

kj  min  (-jjj-,1) 

min  (  .,1) 

(e) 

k. 

J 

l+l 

1 

i+l 

Lemma  1.7  Let  the  notation  b,  B,  p,  and  P  for  binomial  and  Poisson  probabilities  be  as 
in  Chapter  I,  Section  G.  Let  Rj,  R,  d,  dj,  w,  wj,  and  L  be  as  in  Chapter  I,  Section  A; 
assume  that  they  obey  the  restrictions  in  Chapter  I,  Section  B.  As  in  Chapter  I,  Section  G, 
let  p  =  dw/L,  and  define 

(  ^ 


to  =  d 


R-l 

i 


pl(l-p) 


R-l-i 


x=0,...,R-l 


and 


VJt=  dJe  1=0,1, 2,..., 

where,  for  the  purposes  of  this  lemma,  the  aj  are  some  positive  constants.  (The  index  j  is 
redundant  here.)  Then  for  integers  ti  and  12  such  that  0  <  11  <  12: 

l2 

a)  £  «  =  d[®(i2;  R-l,p)  -  B(i,  -1  ;R-1  ,p)]. 


A-10 


b) 


c) 


d) 


1*2 

T  —co  =  -^-[B(i2+l;R,p)  -  ©(ijiR.p)], 
,  ,  i  i  wR  *■ 


i=ij 1+1 


J  ¥jt  =  djtPa^cCj)  -  PCij-Uctj)  ], 
l=M 

h  i  d= 

Zj  —  V.  =  —  [P  (l2+1:  aj)  -P(i,;ou  ]. 

l=lj  l+l  Jl  “j  1  J 


Proof:  It  is  evident  that  cot  =  db(i ;  R-l,p)  and  \|/jt  =  djp(i;aj)  so  formulas  a) 
and  c)  follow  immediately  from  the  definitions  of  B(m;n,p)  and  P(m;  |i)  in  Chapter  I, 
Section  G.  Note  also  that  for  0  <  i  <  R  -  1 


1 


t+1  i 


0)  =  d- 


(R-l)  ! 


(R-l-i)  !  l  !  (l+l) 


Pl  (1— p) 


R-l-i 


Multiplying  and  dividing  the  right-hand  side  by  Rp  and  simplifying  yields 

1  d 

Wv  =  Rp  b(l+1;  R’p) 


l+l 


and  thus 


i2  +  l 


i=R,  i+ 


-  w  =  X  b(m;R,p) 

l+l  1  ^  ,  Rp 


m=ij  +1 


Rp 


[B(i2+l;R,p)  -  B(ij  ;R,p)]. 


Substituting  wd/L  for  p  yields  expression  b).  If  12  ^  R,  B  (i2+l;R,p)  remains  1,  and  thus 
expression  b)  holds  for  all  values  of  li  and  12. 

Similarly 


d. 


-Vji  =  de  Ja-  /(i+l)!  =  of  p(i-+l;ct )  . 

•1  J  J  j  J 

Performing  the  indicated  sum  clearly  yields  expression  d). 


1+] 
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B  .  PROOFS  OF  THEOREMS  IN  CHAPTER  I,  SECTION  G 


This  section  presents  outlines  for  proofs  of  the  theorems  stated  in  Chapter  I, 
Section  G.  The  main  conditioning  arguments  have  been  derived  in  Chapter  I,  Section  E. 
Here,  we  briefly  restate  them  as  appropriate;  see  Section  E  for  justifications.  Actual 
algebraic  simplications  are  not  performed,  but  references  to  lemmas  are  made  where 
appropriate.  Unless  specifically  redefined,  all  symbols  for  parameters,  events,  random 
variables,  and  so  forth,  are  as  in  Chapter  I. 

Homogeneous  (n=l)  and  heterogeneous  (n>  1 )  cases  are  treated  together:  the 
derivations  are  very  similar.  In  the  homogeneous  case,  suppress  the  subscripts  j  on  the 
parameters,  and  consider  cot  as  being  zero  for  i>R. 

In  the  heterogeneous  case  where  some  ccj  are  zero,  some  results  have  a  form 
different  from  the  general  equations.  In  the  proofs  of  the  theorems  we  assume  aj  >0;  in 
Section  4,  below,  we  discuss  the  case  aj=0. 

Independence  conditions,  Poisson  approximations,  and  edge  effects  are  not 
mentioned  in  this  appendix;  they  affect  the  final  results  in  the  ways  discussed  in  Chapter  I, 
Sections  D,  E,  and  F. 

For  the  rest  of  this  appendix,  "Section"  refers  to  the  indicated  section  in  Chapter  I. 

1.  Theorems  1,1  and  1.2:  Expected  Numbers  of  Penetrators  of  Type  j 
Killed 

The  derivation  starts  by  evaluating  the  probability  P(K)  (that  a  specific  barrier 
element  kills  a  specific  penetrator  which  it  can  potentially  detect),  using  the  formula  from 
Section  E.  1  (page  19) 

P(K)  =  £  £  P(KIF  )  P(FhIDAl)P(DAl)  . 

h=l  i=0 

The  P(Fh  IDAO  are  given  by  Lemma  1.5,  P(K  I  Fh)  is  simply  l-(l-kj)^,  and  the  P(DAt) 
are  the  aq  or  (approximately)  vj/jt  defined  in  Section  G.l.  To  evaluate  P(K),  first 
interchange  the  order  of  summation.  The  sum  on  h 
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H 

X  P(KIFh)P(FhIDAv) 
h=l 

results  in  the  P(K  I  DAO  given  in  Lemma  1.6.  Then 

oo 

P(K)  =  y  P(KIDA  )  P(DA  ) . 

i  i 

t=0 

The  results  of  Lemma  1.7  can  be  used  to  evaluate  this  sum  for  attack  protocols  (a),  (d),  and 

kjH 

(e).  In  (d),  the  sum  must  be  broken  into  two  parts:  where  P(KIDAi)  is  kj  or  — j-  , 
respectively.  The  sum  for  protocol  (b)  does  not  appear  to  have  a  closed  form.  In 

H 

evaluating  the  sum  for  (c),  we  note  that  for  i>H,  the  integer  and  fractional  parts  of  — j- 
H 

become  zero  and  — p  and  thus  the  P(KIDAt)  of  Table  A-l  becomes 

IfU 

P(KlDA  )  =  - 

The  tail  sum 

£  £fP<DA) 

i=H 

can  be  evaluated  using  Lemma  1.7  (noting  that  in  the  homogeneous  case  P(DAx)  =  0  for 
t>R,  and  in  the  heterogeneous  case,  letting  equal  1  for  any  |i).  The  rest  of  the  sum 

H-l 

Y  P(KIDA  )P(DA  ) 

l=0  1  1 

for  protocol  (c)  does  not  seem  to  have  a  closed  form. 

The  resultant  P(K)  for  protocols  (b),  (c),  and  (e)  are  stated  in  Theorems  1.1  and  1.2; 
for  protocol  (a),  P(K)  =  djkj.  To  derive  the  overall  E(AR.),  we  First  note  that 

E(£Rj )  =  y  P(S) 
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(Reference  [12]),  and  also 


l-P(S) 

L 

J  PCS  ICX)  -  dx 
o  L 

B 

f]  P(Sb  ICX )  if  the  independence  condition  holds 

b=l 

Nj  (x) 

[l-P(K) ]  . 

Using  Lemma  1.3,  the  integration  on  x  yields 

I  I.+  l 

P(S)  -  [1-P(K)]J  (1-fj)  +  [l-P(K)  ] J 

Ij 

=  [  1— P  (K)]  [l-fjP(K)] 

and  the  resultant  formula  for  E(ARj)  is  as  displayed  in  Theorems  1.1  and  1.2. 

2.  Theorems  1.3  and  1.4:  The  Expected  Number  of  Shots  Fired  at  a  Given 
Penetrator 

First,  consider  one  specific  penetrator  (of  type  j)  and  find  the  expected  value  of  the 
number  of  shots  Ti  fired  at  the  penetrator  by  a  given  barrier  element  that  can  potentially 
detect  it.  By  arguments  in  Section  E  l, 

E(h)  =  X  £  hP(F„IDA  )P(DA  )  . 
h=l  i=0 

The  PfFIDAO  are  given  by  Lemma  1.5;  the  P(DAl)  are  or  \j/jb  for  the  homogeneous  and 
heterogeneous  cases,  respectively.  Interchange  the  order  of  summation;  the  intermediate 
sums 

H 

E(hlDA)=£  hP(FhIDAt) 

1  h=l 

are  derived  in  Lemma  1.6.  Examining  these  results,  it  is  evident  that  for  all  attack  protocols 
the  sums 

oo 

E(h)=X  EfhlDAjPfDAJ 

i=0 

can  be  evaluated  using  Lemma  1.7.  The  results  are  givm  in  Table  A-2. 


P(S)  = 
P(S)  = 

P(SICX)  = 
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Table  A-2.  EXPRESSIONS  FOR  E(h) 


Attack  Protocol 

Homogeneous  Case 

Heterogeneous  Case 

(a) 

d 

dj 

(b) 

wR  1  V  L  '  J 

Hd:  -aj 

— -  (1— e  J) 

(c) 

“  [i-(i-"d)R] 
wR  L 

Hdj  _a. 

(1— e  J) 
aj 

(d) 

If  H>R,  E(h)=d 

If  H<R-1 

E(h)  =  dB(H-l;R-l,p) 

HI 

+— -[l-B(H;R,p)  ] 
wR 

where  p=wd/L. 

d}nn-v,a}) 

Hdj 

+  ttj  [l-P(H;aj )  ] 

(e) 

wR  [1-”-  WL)R] 

d  -<x 

—  (1-e  J) 

aj 
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To  compute  E(y),  the  expected  number  of  shots  fired  at  a  given  penetrator,  p,  we 
note  that 


B 

E(y)  =  X  E(yb)  . 

b=l 

where  is  the  expected  number  of  shots  fired  at  the  penetrator  by  barrier  element  b.  If  b 
can  potentially  detect  p,  y^  has  the  same  distribution  as  H;  if  b  cannot  potentially  detect  p, 
y^  is  zero  with  certainty.  Given  that  p  crosses  at  point  x,  Nj(x)  barrier  elements  can 
potentially  detect  p.  Thus 

B 

E(ylQ)  =  SE(ybiq) 

b=l 


=  Nj  (x)  E(h)  +  (B-Nj  (x) )  •  0 

Integrating  on  x  and  using  Lemma  1.3, 

L 


E(y)  = 


0 


E(ylCx)j-  dx 

L 


=  (l-fj)IjE(fi)+  fj  (Ij+DECfi)  . 

Simplifying  and  recalling  that  (by  definition)  Ij  +fj  =  wjB/L, 

w.B 

E(y)  =  -j~-  E(h). 

In  the  homogeneous  case,  the  wj  and  L  will  cancel  terms  from  E(K).  The  final  results  for 
E(y)  are  as  stated  in  Theorems  1.3  and  1.4. 


3.  Theorems  1.5  and  1.6:  The  Probability  Distribution  of  the  Number  of 
Shots  Fired  at  a  Given  Penetrator 

The  proof  of  Theorems  1. 5a  and  1. 6a  (when  attack  protocol  (a)  is  used)  has 
essentially  been  given  in  Section  E.2.  We  note  here  that  Lemma  1.3  is  used  in  integrating 
over  x.  The  algebraic  result  for  P(y=y)  in  that  section  is  the  same  as  the  result  stated  in 

Theorem  1. 6a.  Note  that  if  f  =0,  P(y=0)  =  1-fj  dj  and  P(y)  =  f j  d j .  This  makes  sense:  a 

penetrator  receives  one  shot  if  and  only  if  it  is  detectable  by  some  attacker  (which  occurs 
with  probability  fj)  and  that  attacker  then  actually  detects  the  penetrator. 
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For  the  other  attack  protocols,  we  assume  that  wj  <  L/B  for  all  j=l,...,n.  Lemma 
1.4  reduces  the  problem  to  finding  the  probability  P(Fy)  that  a  specific  barrier  element  (3 
fires  exactly  y  shots  at  a  specific  penetrator  p  that  the  barrier  element  can  potentially  detect. 
The  quantity  y  can  range  from  zero  through  H,  the  input  number  of  shots  per  shooter.  In 
attack  protocol  (e),  set  H=l.  The  formula  in  Section  E.l 

oo 

p(Fh )  =  P(FhID)P(D)  +  X  P(Fh!DA  )P(DA  ) 

i=0  1  1 

yields 

P (F0 )  =  l^dj  +  £  P(F0IDA  )P(DA  ) 

1=0 

and 

P(FV)  =  £  P (F  i  DA  )  P(DA  ),  y=l,...,H  , 

y  >  t  i 

1=0 


using  y  rather  than  h  as  an  index.  The  P(Fy  IDAO  are  given  in  Lemma  1.5;  the  P(DAt)  are 
the  on  or  \j/jl  defined  in  Section  G.l.  Finding  the  P(Fy)  then  becomes  a  matter  of 
evaluating  the  indicated  sums  above.  Looking  at  the  Lemma  1.5  results,  we  see  that  the 
results  of  Lemma  1.7  can  be  applied  in  a  straightforward  manner  to  attack  protocols  (d)  and 
(e)  for  the  case  y=l.  Since  in  these  protocols,  a  given  barrier  element  will  never  fire  more 
than  one  shot  at  a  given  penetrator,  P(Fo)  =  l-P(Fj).  The  results  are  as  stated  in  Theorems 
1.5,  and  1.6  d  and  e.  The  formulas  for  attack  protocol  (b)  seem  to  be  impossible  to 
simplify. 

The  interesting  case  is  attack  protocol  (c).  From  now  on  we  use  h  rather  than  y  as 
the  number  of  shots  fired.  From  Lemma  1.5, 


0  for  all  other  h. 
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To  evaluate 


£  P(F^IDA)P(DA  ) 
1=0  1  1 


for  given  h,  we  must  determine  which  i  yield  nonzero  P(FhlDAi).  We  define  the  following 
two  series  of  sets,  indexed  by  h=0,...,H,  where  l  can  assume  only  nonnegative  integer 
values.  Let 


L](h) 
L2  (h) 


1 1 1  h  - 

1  All 

Lt+i  JJ 

1  H  1 

j  1 1  h  = 

Li+iJ 

For  a  fixed  nonnegative  integer  h,  h 


(i.e.,  i  e  Lj  (h) )  if  and  only  if 


h  < 

Rearranging,  if  h>l,  this  is  equivalent  to 

H 


H 

l+l 


<  h+1 


H 


.  -1  <  i  <  -1 

h+1  h 


Similarly,  it  can  be  shown  that  if  h>2,  then  h 


i.e.,  ig  L^Ch),  if  and  only  if 


Therefore,  for  h>2, 


<  l 


< 


H 

h-1 


-1  . 


L.(h)  -  {'^r " ' s  r-  *} 

The  cases  h=0  and  h=l  must  be  treated  specially.  Note  that  I3J  =0  if  and  only 


if  H<i+1.  Throughout,  we  are  assuming  that  H  is  a  positive  integer.  Thus 
if  and  only  if  t>H,  and 


=0 
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L,(0)  =  {  H,  H+l.  H+2,  ...} 

By  inspection,  1.2(0)  is  the  empty  set,  and  L2(1)=Li(G);  Lj(i)  can  be  found  using  the 
general  formula. 

It  is  evident  that  for  any  given  h,  Li(h)  and  L2O1)  are  disjoint,  and  that  for  hi*h2, 
L i (hi)  and  L](h2)  are  disjoint,  as  are  L2(hi)  and  L2(h2).  Also,  L2(h+1)  =  Li(h)  for 
h=0,...,H-l.  Some  of  the  Li(h)  and  L2(h)  may  be  empty. 

We  can  find  integer  bounds  on  the  i.  Let  k  be  any  integer  and  x  any  real  number 
(possibly  integer).  It  is  clear  that  k<x  if  and  only  if  k  <  LxJ,  and  it  can  be  easily  shown 

(treating  the  cases  for  integer  and  noninteger  x  separately)  that 

x  <  k  if  and  only  if  LxJ  +  1  <  k  . 

Therefore  the  sets  Lj(h)  become 

L , (0)  =  {  H,  H+l,  H+2  ....} 

L2(0)  =  <t> 

l,  (i)  =  { 1 1  LyJ  <  i  <  H-i} 

L2(l)  =  {  H,  H+l,  H+2,  ...} 

and  if  H>2, 
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Thus  the  sum 

oo 

X  P(FhIDA  )  P(DA  ) 

i=0  1 

used  in  the  computation  of  the  P(Fh)  becomes 

X  (h+1  -  — )P(DA  )  +  X  (—  -  [h— 1  ] )  P(DA  )  . 

leLj(h)  l+i  1  ieL2(h)  i+l 

Since  P(DAl)  is  or  \|/ji,  we  can  evaluate  this  sum  using  Lemma  1.7,  with  the  appropriate 
integer  limits  for  ILj  (h)  and  L2(h)  used  as  limits  of  summation.  The  set  Lj(h)  is  empty  if 
the  corresponding  expression  in  Lemma  1.7  is  negative,  so  we  use  the  positive  parts  of 
such  expressions.  The  formulas  in  Theorems  1. 5c  and  1. 6c  follow  immediately. 

4 .  The  Heterogeneous  Case  Where  aj=0 

When  the  Poisson  approximation  to  the  binomial  is  used  to  determine  the  probabili¬ 
ties  P(DAi),  ctj,  which  is  defined  by 
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If  aj=0,  P(DAo)=dj  (which  equals  P(D) )  and  P(DAt)=0  for  all  l>l.  The  general  formulas 
of  Chapter  I  still  hold;  the  appropriate  P(FhlDAt)  can  be  substituted  in  and  the  results 

follow  immediately. 
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PROOFS  OF  RESULTS  IN  CHAPTER  II 


Appendix  B  is  organized  similarly  to  Appendix  A.  The  first  part  proves  three 
lemmas  used  in  the  derivations.  Lemmas  II.  1  and  II. 2,  stated  in  Chapter  II,  Section  C, 
concern  the  distribution  of  the  events  Wm.  Lemma  II. 3  derives  formulas  for  certain  sums 
and  weighted  sums  of  the  0jt  (i.e.,  P(DA^)  ).  The  second  part  of  the  appendix  provides, 
for  each  theorem  of  Chapter  II,  Section  E,  proof  details  to  complete  the  derivation  steps 
outlined  in  Chapter  II. 

Unless  specifically  redefined,  notation  for  all  events,  random  variables,  parameters, 
and  so  forth,  is  identical  to  that  used  in  Chapter  II. 


A.  THE  LEMMAS 


Let  Ij  be  the  integer  pan  and  fj  the  fractional  part  of  wjB/L.  Let  the 


events  V  and  Wm  be  as  defined  in  Chapter  II.  Then,  ignoring  edge  effects. 


P(Wj  IV) 

j 

P(WmIV) 

J 


h  q-*j> 

(Wj  B/L) 
(Wj  B/L) 


P(Wm  I  V)  =  0  for  all  other  m. 

Proof:  Let  the  specific  barrier  element  under  consideration  be  denoted  as  barrier 
element  P  and  the  specific  penetrator  (of  type  j)  as  penetrator  p.  Note  that  p  is  an  integer 
between  one  and  B.  Barrier  element  P's  position  on  the  barrier  line  is,  by  Assumption  1) 


x3  =  •§-(|3-'f)  • 


Let  x  be  the  crossing  point  of  penetrator  p.  By  Assumption  4),  if  V  occurs,  x  must  lie 
within  wj/2  of  xp.  We  ignore  edge  effects  by  assuming  that  the  distances  of  xp  from  the 
edges  of  the  barrier  (points  0  and  L)  are  at  least  Wj/2.  (If  Wj  <  L/B,  this  can  easily  be 
demonstrated.)  Therefore,  the  interval  of  width  wj  centered  on  point  xp  does  not  overlap 
the  edge  of  the  barrier.  The  conditional  distribution  of  x  given  V  is  therefore  uniform  on 
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[xp-wj/2,  xp+wj/2  ].  We  are  interested  in  the  total  number  of  barrier  elements  (including 
(3)  that  p  can  potentially  detect,  given  V.  This  depends  on  the  specific  value  x  of  p's 
crossing  point.  Let  Nj(x)  be  this  number  of  barrier  elements.  Nj(x)  is  the  cardinality  of  the 
set 


{ b  1 1  xfe  -x  I  <  wj  /2 }  , 


where 


is,  by  Assumption  1),  the  location  of  barrier  element  b.  In  reality  b  can  asume  only  integer 
values  between  one  and  B,  inclusive,  but  we  ignore  edge  effects  by  letting  b  assume  any 
integer  value.  For  any  given  x,  the  condition 

I  xb  -x  1  <  wj/2 


is  algebraically  equivalent  to 


/ 

(  „  y 

X 

to  s 
r 

v _ 

<  b  ^  2  B 

X  W- 

+  B 

Let  us  denote  the  left  endpont  of  this  interval  by  b()(x);  the  right  endpoint  is  then  bo(x)  + 


(wjB/L).  If  bo(x)  is  not  an  integer  (and  given  that  %  is  uniform,  bQ(>0  will  be  integer 
with  probability  zero),  then  Nj(x),  the  number  of  integer  points  b  in  the  interval 

b()(X)’  b0(x)  +  wjB//L 

can  be  shown  by  the  same  methods  used  in  the  proof  of  Lemma  1.3  to  be 

Nj  (x)  =  Ij  +  [<b0  (x)  >  +  fjj  # 


It  is  clear  that 


[_<bQ  (x)  >  +  fj J  - 


if  <±>q  (x)>  <  1— fj 


1 


if  <bQ  (x)> 


>  1-fj. 


Thus,  the  random  variable  Nj(x)  will  assume  only  the  values  Ij  or  Ij  +1,  and  will  equal  Ij+1 
if  and  only  if  <bQ(x)  >  >  1-fi  .  (Wm  occurs  if  and  only  if  NjfX)  =  m.)  We  derive  the 
probability  of  this  event. 


Given  V,  the  crossing  point  x  is  uniformly  distributed  on 


xp-wj/2’  xp+wJ/2 


1 


w 


Thus,  b^(x)  =  —  +  B( - -)  is  uniformly  distributed  on 

0  2  L  2L 


1  B  .  1  B 

’TV1'  2+  L  Xpj 


L 

Substituting  in  g"  \P  _  2  /  for  xp  and  simplifying,  this  interval  becomes 


Recall  that  (3  is  a  positive  integer.  The  above  interval  is  the  disjoint  union  of  the  sequence 
of  intervals 

[P  -  (wjB/L),  [3  -  Ij] 

(P-Ij.p-dj-l)] 

(p-(lj-l).p-(Ij-2)] 

(P-l.P]  . 

The  first  interval  has  length  fj,  the  others  have  length  unity. 

Let  the  random  variable  Z  be  defined  as  the  fractional  part  of  bgCx)  ;  i.e., 

2  =  <bQ(x)>  . 
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We  wish  to  determine 


P(2">  1  -  fj). 


If  bgCx)  lies  in  one  of  the  Ij  unit  intervals  (which  happens  with  probability  Ij/(Ij  +  fj ) ),  z 
will  be  uniformly  distributed  on  [0,1]  and 


P(z  >1  -  fj)  =  fj  . 

If  b^(x)  lies  in  the  interval  [(3-Bwj/L,  p— Ij]  (which  happens  with  probability  fj/(Ij  +  fj) ), 
then  z,  the  fractional  part  of  bg(x) .  will  certainly  be  at  least  1  — fj.  Therefore, 


P(Z  >  1-f, )  = 


!j+fi 


fj 


h«i 


fj(Ij+D  fjOj+1) 
Ij  +  fj  "  (WjB/L) 


By  our  previous  reasoning,  this  is  the  probability  that  Nj  (x)  =  L  +1 ,  and  thus 

f.di+1)  Ij  (1-9 


P(Nj  (50 


j '  J 


Ij)  1  (WjB/L)  (WjB/L)  ’ 


which  is  the  content  of  the  lemma. 


Lemma  jl, 2.  Let  Ij,  fj,  and  Wm  be  as  defined  in  Chapter  II  (and  in  Lemma  II. 1). 
Then,  ignoring  edge  effects: 

P(W,  )  =  1-fj 
PCWi  +  i  )  =  ^ 

j 

P(Wm)  =  0  for  all  other  m. 

Proof.  Again,  call  the  specific  penetrator  under  consideration  penetrator  p, 
assumed  to  be  of  type  j.  By  Assumption  1),  barrier  element  b  is  positioned  at 
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and  by  Assumption  2),  p's  crossing  point  5T  is  uniformly  distributed  on  [0,L],  By 
Assumption  4),  for  any  given  x,  Nj(x),  the  number  of  barrier  elements  that  p  can 
potentially  detect,  is  the  cardinality  of  the  set 

w. 

{b  I  I  xb-x  I  <  -j- }  . 

We  ignore  edge  effects  by  considering  any  integer  b,  even  though  b  can,  in  fact,  only  be 
between  1  and  B.  If  wj  <  L/B  it  is  easily  shown  that  for  any  x  e  [0,L],  the  above  set  will 
contain  no  integer  points  b  outside  of  [  1  ,B] .  Exactly  as  in  the  proof  of  Lemma  II.  1,  the 
condition  I  xb  -x  I  is  algebraically  equivalent  to 

Bw 

bQ(x)  <  b  <  bQ(x)  +  , 


where 


1  B  wj 
b0<x)  =  2  +  L(X  "  —  }  ■ 

Also,  as  in  Lemma  II.  1,  the  number  of  integer  points  b  in  [bo(x),  bQ(x)  +  Bwj  /L]  is  Ij+1 


if  and  only  if  <bO(x)>  >  1  -  fj.  Now,  however,  x  is  uniformly  distributed  on  [0,L]  (rather 
than  [xp  -  wj  /2,  xp  +  wj  /  2]  )  and  thus  bglx)  is  uniformly  distributed  on 


Bw,  i  Bw- 

0  -  V  >  •  2  (1  -  TT*  +  B1  • 


Since  B  is  a  positive  integer,  this  interval  is  composed  of  B  intervals  of  length  1.  As  bp  (x) 
ranges  over  this  interval,  the  fractional  part  of  bq  (x)  assumes  in  turn  every  value  between  0 
and  1,  taking  on  each  value  B  times.  Thus,  <bg  (x)  >  is  uniformly  distributed  on  [0,11  and 
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P(<bQ(x)>  >  1-fj )  =  f.  . 

But  by  previous  reasoning,  this  is  the  probability  that  p  can  potentially  detect  Ij  +1  barrier 
elements. 


Lemma  II. 3.  Let  Ij  and  fj  be  as  in  Lemmas  II.  1  and  II. 2.  As  in  Chapter  II, 
Section  E,  let  the  probabilities  PCDAv)  be  denoted  by  0ji,  where 

L-l-i 

T  ■  n-f.t  A.  I  ‘  I  H-A 
B 

I  .-l 

n.  +n  f.  a.  I  1J  Ih1  n-A. ) 

B 

and 


0ji  =  -g-,j(1-fj)di(Iji"1)  dj  oV' 

4<Ij+,)fJd)('')d'J0^J)'rl  for  l=0,...,lj-l 


1+  1 

0  =  J_  (Ij+l)fj  d  J  . 

Jlj  B  J  J  J 


Then 


a)  For  0<li<t2SIj-l 


I,(l-f,)d, 


Vn  =  J.: 


_  0 

t=l.  Jl 


[E(i2;  Ij  —  1 ,  d j )  -  Bdj-ljlj-l.dj)] 
(I  +Df  d 

+  B  [B(i2;  Ij,  dj)  -  B(i j-LIj.dj)] 


b)  For  0  <  ii  <  12  <  Ij—  1 


£  ITT  9j>  =  -g- O-fjjf - Bfi,;r,dj  >] 
1 

+  i-f  [ B(i  +1;I  +l,d  )-B(i  ;I  +l,d  )]. 

J  J  J  1  J  J 
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c)  For  0  <  l!  <  Ij 


I,  (1-fj )  d 

Xe  =  J  RJ  J  -1;  I- -1,  d- )  ] 

i=ij  ji  a  3 

(I;+l)  fjd: 

+  J  -  J  J  [1-B(i,-1;  Ij ,  dj)  ] 


d)  ForO  <  ii  <  Ij 


-J-0. 

l+l  J>- 


♦£fj[1-*V,Jtl-dJ)1- 


where  the  notation  for  binomial  probabilities  is  as  in  Chapter  I,  Section  G. 

Proof.  For  simplicity  in  notation  we  suppress  the  subscript  j  on  dj,  fj,  and  Ij.  Let 

*1  =  ^-1  (1-0  d 


and 

a2  =  X(i+i)fd  . 

Note  that  for  t  =  0,...,I-1 

6  =  a,  b(i  ;I-1,  d)  +  a7ib(i;I,d) 

ji  1  z 

and 

0j,  =  aj  1>  (I  ;  I,  d)  . 
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Also,  for  i  = 


7+T  b(l;H,d)  -  ^ 


x 

Id 


I ! 


(I— [i  +1]  )  !  (i+l)  ! 


dl+l(l-d) I_[l+1] 


=  tv  b(i+l;I,d) 
Id 

Thus,  — j-  b(i;I-l,d)  =  -^-(1-0  b(i+l;  I,d)  . 


For  i  =  0,...,I 


i+l 


b  (i;I,d)  =  -*■ 


I!  dl 


l+l  (I-i)  !  i  ! 

1  (I+l)  ! 


(I+l)d  ([I+l]-[t+l]  )  !  (i+l)  ! 
1 


-dl+W~l 


(I+l)d 


b(i+l;I+l,d) 


Thus,  -i-  a,,  b(i;I,d)  =  A-fb  (i+l;  I+l,d)  . 
1+ 1  ~  B 


Thus,  for  0  <  li  <  12  <  Ij-1 
l2  ’’2 

£  0lt  =  £[a1b(i;I-l,d)  +  a2b(i;I,djl 


ii  '  l 

l=li  l=li 


and 


X  T+t6!!  =  X  +d-flb(i+l;I,d)  +  X  X<i+1;  T+l,d). 

i=i j  1  J  t=i|  °  i=ij  D 


Formulas  a)  and  b)  follow  forthwith  from  the  definition  of  the  cumulative  binomial. 
Assume  that  \\  <1-1.  Then 

I  I-i  I 

X  e„  =  X  ai  b(l;  I-l.d)  +  X  ao  b(l’ 

1=1 !  l=ll  l=ll 
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and 


Formulas  c)  and  d)  follow  forthwith. 

If  ii  =  I,  tlrn 

I 

X  e,  =  a,  e|T  =  a0b(I;I,d) 

1=1  *-  JA  ^ 

1  ll 

If  ii  =  I  is  substituted  in  formula  c),  the  First  term  becomes  zero,  and  the  second  term  is 

a2  [l-B(I-l;I,d)]  =  a2b(I;I,d)  . 

In  formula  d),  substituting  I  for  ij  yields 

0  +  4-f[l-B(I;I+l,d)], 

r> 

which  equals 

4-  fb(I+l;I+l,d), 

D 

or 


Thus  formulas  c)  and  d)  are  true  for  all  i,  between  0  and  I,  inclusive. 

B  .  PROOF  DETAILS  OF  THEOREMS 

The  discussion  in  Section  C  of  Chapter  II  has  presented  the  overall  steps  in  the 
derivations  of  the  desired  quantities.  First,  quantities  qj,  ej,  and  <j>hj  are  derived  for  one 
penetrator  and  one  barrier  element.  These  are  then  combined  as  explained  in  Chapter  II, 
Section  C.2.  As  derived  in  Section  C.l,  the  formulas  for  qj,  ej,  and  <J>hj  are 
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[ :  =  t  ^  P(KlF^)  P(FhIDAt)  PHDA^ 


h=0  i=0 
=  gP(KlDAt)P(DAv) 

1  Ij 

=  £  ^  hPCF^DAOPCDAO 
h=0  i=0 

=  ^E(KDAi)P(DAt) 

i=0 


and 


0Oj  =  1 


wd- 
j  j 


i, 

t 


i=0 


PCFolDA)1'  .DA) 


p(^i|DAi)P(DA  ),  h  = 

i=0 


The  P(FhlDAt)  are  given  by  Lemma  1.5  in  Appendix  A,  the  P(DAt)  are  simply  the  0jL  of 

Lemma  II. 3,  and  the  P(KIDAi)  and  EChlDA^)  are  derived  in  Lemma  1.6.  Examining  these 

quantities,  we  see  that  most  of  the  P(KIDAt),  E(H IDAO,  and  P(FhlDAt)  are  either  constants 
with  respect  to  i  or  involve  a  denominator  of  i+l.  Lemma  II. 3  can  then  be  used  to  evaluate 
the  above  indicated  sums.  In  attack  protocols  (b)  and  (c),  we  have  not  been  able  to  derive 
closed  forms  for  qj  and  4>hj»  except  as  noted  in  Theorems  I1.3b  and  c. 


H  H  = 

min  ( — V,  1),  and  P(KlDAi)  =  k:min(— i-,  1)  .  To  evaluate 


The  interesting  case  is  attack  protocol  (d),  where  PfFilDAt)  =  E(fnDAt)  = 

L 

i+i  ’ 

i. 

e.  =  { 


u+1 


min 


i=0 


H. 

J 


A 


w-i’1. 


e. 

ji 
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we  first  note  that  the  minimum  is  less  than  1  if  and  only  if  t>Hj.  If  Hj>Ij+l,  for  all  values 
of  i  with  positive  0^,  the  EfhlDAO  term  will  equal  1.  Lemma  II. 3  formula  c)  then  yields 
(treating  B(-l;n,p)  =  0) 


w.d. 
_  J  J 

L 

If  Hj<Ij,  the  formula  for  ej  becomes 


These  sums  can  be  evaluated  by  formulas  a)  and  d)  of  Lemma  II. 3,  respectively.  When 


multiplied  by  kj,  ej  becomes  qj,  which  is  presented  in  the  statement  of  Theorem  II.  Id. 
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